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ABSTRACT 

We present Chandra monitoring data for six gravitationally lensed 
quasars: QJ 0158-4325, HE 0435-1223, HE 1104-1805, SDSS 0924+0219, 
SDSS 1004+4112, and Q 2237+0305. We detect X-ray microlensing variability 
in all six lenses with high confidence. We detect energy dependent microlensing 
in HE 0435-1223, SDSS 1004+4112, SDSS 0924+0219 and Q 2237+0305. We 
present a detailed spectral analysis for each lens, and find that simple power-law 
models plus Gaussian emission lines give good fits to the spectra. We detect in- 
trinsic spectral variability in two epochs of Q 2237+0305. We detect differential 
absorption between images in four lenses. We also detect the FeKa emission 
line in all six lenses, and the Ni XXVII Ka line in two images of Q 2237+0305. 
The rest frame equivalent widths of the FeKa lines are measured to be 0.4-1.2 
keV, significantly higher than those measured in typical active galactic nuclei of 
similar X-ray luminosities. This suggests that the Fe Ka emission region is more 
compact or centrally concentrated than the continuum emission region. 

Subject headings: accretion, accretion disks — black hole physics — gravitational 
lensing — quasars: individual (QJ 0158-4325, HE 0435-1223, HE 1104-1805, 
SDSS 0924+0219, SDSS 1004+4112, Q 2237+0305) 



1. Introduction 

The concept of quasar microlensing, the effect of gravitational lensing by stellar mass 
objects in a lens galaxy on background quasars, was proposed almost immediately (e.g., 
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Chang & Refsdal 1979; Gott 1981) after the discovery of the first gravitational lens (Walsh 
et al. 1979). However, it was only ten years later that the effect was observationally detected 
in Q 0957+561 (Vanderriest et al. 1989) and Q 2237+0305 (Irwin et al. 1989). Recently, 
quasar microlensing has become an important tool for probing the structure of quasar ac- 
cretion disks and distinguishing models for their non-thermal emission (see the reviews by 
Wambsganss 2006; Kochanek et al. 2007). The ensemble of microlensing stars produces a 
magnification pattern in the source plane, and as the source moves across this pattern due 
to relative motions between the source, the lens, the stars in the lens, and the observer, the 
source magnification varies with time. More importantly for studying accretion disk models, 
the (time varying) microlensing magnification depends on the size and spatial structure of 
the emission source in the central engine of quasar. Larger emission regions smooth the 
magnification patterns more heavily and will have less microlensing variability than smaller 
emission regions. Additionally, the magnification diverges on caustics, which makes it pos- 
sible to resolve emission regions arbitrarily smaller than the Einstein radius of the stars. 
Unlike almost any other astrophysical probe, microlensing becomes a more powerful tool as 
the source becomes smaller. As a result, quasar microlensing is a unique tool for probing 
the innermost structure of accretion disks. 

There are two major challenges in utilizing this approach: first, to measure the mi- 
crolensing light curves and second, to interpret the curves. There are several effects con- 
tributing to the light curves of quasar images and their variability. First, the average fluxes 
of the images are affected by the macro lens model, milli-lensing by substructures in the 
lens (e.g., Kochanek & Dalai 2004; Cohn & Kochanek 2004; Keeton & Moustakas 2009), 
extinction by dust (e.g., Falco et al. 1999; Eliasdottir et al. 2006), and X-ray absorption 
by gas in the lens (e.g., Dai et al. 2006; Dai & Kochanek 2009). Second, they vary with 
time due to a combination of microlensing and the intrinsic variability of quasars. Of these 
effects, the intrinsic variability is usually relatively easy to remove by measuring the time 
delays between the images and then shifting the light curves. Many of other effects are 
also energy/wavelength dependent, but the best way to distinguish microlensing variability 
from other effects is through long term monitoring programs at the relevant energies. This 
approach also has the advantages of ameliorating systematic uncertainties due to the mean 
stellar mass in the lens galaxy. Limited by observational resources, early microlensing studies 
were confined to either single epoch, multiple wavelength (e.g., Falco et al. 1996; Mediavilla 
et al. 1998; Agol et al. 2000), or multiple epoch, single wavelength studies (e.g., Wozniak 
et al. 2000). Only recently have there been intensive studies spanning both the time and 
wavelength (e.g., Poindexter et al. 2008; Chen et al. 2011; Blackburne et al. 2011; Mosquera 
et al. 2011; Muhoz et al. 2011). 

Theoretically, the importance of the clustering of caustics on interpretation of quasar 
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microlensing light curves was realized quite early (e.g., Wambsganss et al. 1990). Because 
of the high optical depth of quasar microlensing in lens systems, quasar microlensing light 
curves are generally affected by many stars, thus requiring theoretical models constructed 
for these complicated conditions. The only real solution has been to simply model the data 
using variants of the Bayesian Monte Carlo method (Kochanek 2004). This technique and 
its variants have led to unique constraints on the emission sizes of accretion disks in different 
wavelengths (Morgan et al. 2008; Eigenbrod et al. 2008; Chartas et al. 2009; Dai et al. 2010; 
Mediavilla et al. 2011, Poindexter & Kochanek 2010b; Blackburne et al. 2011), and on the 
stellar/dark matter distribution of the foreground lensing galaxies (e.g., Mediavilla et al. 
2009; Poindexter & Kochanek 2010a; Bate et al. 2011; Pooley et al. 2012). 

We are particularly interested in using microlensing to study the non-thermal, X-ray 
emission regions of quasar accretion disks. The X-ray continuum emission of quasars is 
observed to follow a power-law, and it is believed to be generated through the unsaturated 
inverse Compton scattering of UV disk photons on hot electrons in a corona above the disk 
(see the review by Reynolds & Nowak 2003). The origin of the corona, its geometrical size 
and detailed structure are still topics under debate. In addition to the X-ray continuum, the 
disk-corona model also predicts the existence of a reflection component including emission 
lines (e.g., FeKa) formed through irradiation of the colder accretion disk by the continuum 
(e.g., George & Fabian 1991). However, the location of this reflection component, whether 
it comes from the disk, the disk- wind, or the "torus", has been debated for more than 
two decades. In principle, with enough observations, the quasar microlensing technique can 
resolve all these issues. 

Microlensing of X-rays was first tentatively suggested in RX J0911. 4+0551 (Morgan 
et al. 2001). Because of the lack of multi-epoch observations and the large intrinsic X-ray 
variability in quasars, later studies focused on microlensing of the FeKa line (Chartas et 
al. 2002, 2004; Dai et al. 2003; Ota et al. 2006), where dramatic changes in the FeKa line 
equivalent widths were most easily interpreted. Studies of the X-ray continuum resumed 
with the report of the dramatic flux ratio differences between the X-ray and optical bands 
in RXJ1131— 1231 (Blackburne et al. 2006) and the single epoch, optical-to- X-ray flux ratio 
comparison for a sample of lenses by Pooley et al. (2006). Short Chandra monitoring pro- 
grams of a few lenses start to emerge at this stage (e.g., Chartas et al. 2009), and combining 
them with the Bayesian Monte Carlo microlensing analysis technique (Kochanek 2004), we 
have generally found that the X-ray continuum is emitted on scales of roughly 10 r g (Mor- 
gan et al. 2008; Chartas et al. 2009; Dai et al. 2010; Blackburne et al. 2011), 7-10 times 
smaller than the observed optical emission region (Morgan et al. 2010). However, the cur- 
rent lower limits on the sizes of the X-ray emission regions are not well constrained, because 
strong lower limits require more densely sampled light curves than any presently available. 
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Nonetheless, these microlensing results have ruled out disk-corona models where the corona 
covers extensive areas of the optical accretion disk. The large equivalent widths of the Fe Ka 
line measured in a few lensed quasars suggest that a significant fraction of the FeKa line 
flux originates from a region smaller than the X-ray continuum emission region (Dai et al. 
2003; Popovic et al. 2003, 2006). However, these studies were limited by the lack of deep, 
better sampled X-ray microlensing light curves. 

In this paper, we report results from monitoring six gravitationally lensed quasars with 
Chandra: QJ 0158-4325 (Maza et al. 1995; Morgan et al. 1999), SDSS 0924+0219 (Inada 
et al. 2003a), SDSS 1004+4112 (Inada et al. 2003b), HE 0435-1223 (Wisotzki et al. 2002), 
HE 1104-1805 (Wisotzki et al. 1993), and Q 2237+0305 (Huchra et al. 1985)0. This is 
the largest quasar microlensing X-ray monitoring program to date. QJ 0158—4325 and 
HE 1104—1805 are doubly imaged while the other four systems have four bright images. Five 
of the foreground lenses are galaxies, while SDSS 1004+4112 is lensed by a galaxy cluster. 
Since microlensing in the optical bands has been extensively studied in these systems (e.g., 
Morgan et al. 2006; Kochanek et al. 2006; Lamer et al. 2006; Fohlmeister et al. 2007, 
2008; Eigenbrod et al. 2008; Faure et al. 2009; Poindexter & Kochanek 2010a, 2010b; 
Ricci et al. 2011; Blackburne et al. 2011), our X-ray microlensing light curves provide 
complementary data to constrain non-thermal emission processes in accretion disks. For 
example, a simple comparison of the microlensing variability between the soft and the hard 
X-ray bands in Q 2237+0305 has already led us to conclude that the hard X-ray continuum 
emission region is smaller than that of the soft in this quasar, suggesting a temperature 
gradient in the X-ray emitting corona (Chen et al. 2011). In our original proposals, only 
the observations of Q 2237+0305 and RXJ1131— 1231 (which will be discussed in Chartas 
et al. 2012) were explicitly designed to have the count rates needed to obtain accurate 
measurements in multiple energy bands. 

We present the Chandra observations and the data reduction techniques in Section 2, the 
image analysis in Section 3, and the spectral analysis in Section 4. The microlensing analysis 
is in Section 5, and the discussion and conclusions in Section 6. Throughout this paper, we 
assume a flat ACDM cosmology with = 0.73, Qm — 0.27, and H = 70kms _1 Mpc -1 . 



1 Basic information for each lens system (e.g., the lens/source redshifts, celestial coordinates, and rela- 
tive image positions) can be found on the CfA-Arizona Space Telescope Lens Survey (CASTLES) website 
( http: / / cfa-www.harvard.edu/castles / ) . 
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2. Observations And Data Reduction 

All observations presented in this paper were made with ACIS (Garmire et al. 2003) 
on the Chandra X-ray Observatory (Weisskopf et al. 2002). With an on-axis point spread 
function (PSF) of 0'.'5, Chandra is the only X-ray telescope that can resolve lensed images 
typically separated by 1-2". The majority of the observations are from our large Chandra 
Cycle 11 program to monitor seven gravitational lenses. Here, we report on six of the lenses, 
QJ 0158-4325, SDSS 0924+0219, SDSS 1004+4112, HE 0435-1223, HE 1104-1805, and 
Q 2237+0305. Although the absorption corrected count rates of Q 2237+0305 were published 
in Chen et al. (2011), we present the detailed spectral analysis in this paper. The analysis of 
RXJ1131 — 1231 will be presented in a companion paper (Chartas et al. 2012, in preparation). 
We also analyzed the archival Chandra data for the six lenses. We list the basic information 
for these observations, such as observation dates and exposure times, in Tables [T] — [6j 

We reprocessed all data using the CIAO 4.3 software tools by first removing the pixel 
randomization and applying a sub-pixel algorithm for event positions (Tsunemi et al. 2001). 
We then filtered all events using the standard ASCA grades of 0, 2, 3, 4 and 6 with energies 
between 0.4 and 8.0 keV (observer's frame), and further separated the events into soft and 
hard bands for subsequent analyses. We selected the energy boundary between the soft and 
hard band for each lensing system to balance the soft and hard band count rates (see e.g., 
Table [Q in order to produce the best statistics for detecting differences between the two 
energy bands. For QJ 0158-4325, SDSS 0924+0219, and SDSS 1004+4112 we split the soft 
and hard bands at 3.0 keV (rest frame), for HE 0435-1223 and Q 2237+0305 we split the soft 
and hard bands at 3.5 keV (rest frame), and for HE 1104—1805, we split the soft and hard 
bands at 4.0 keV (rest frame). In Chen et al. (2011), we used (observed frame) full, soft, and 
hard bands at 0.2-8.0 keV, 0.2-2.0 keV, and 2.0-8.0 keV respectively. In the current paper, 
we define the full band as 0.4-8.0 keV, and divide the soft and hard bands at 1.3 keV (3.5 keV 
in rest frame). This energy division provides better statistical results than our definitions in 
Chen et al. (2011), but it alters none of our conclusions. The pile-up effect is negligible for 
all our sources (<2% in all images), since they are generally faint. Nevertheless, we used the 
1/2 subarray mode in our observations of QJ 0158-4325, SDSS 0924+0219, HE 0435-1223, 
HE 1104—1805, and Q 2237+0305 in order to minimize any residual pile- up effects. In our 
program, RXJ1131— 1231 is the only case where the pile-up effect needs to be extensively 
modeled (Chartas et al. 2009). 
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3. Imaging Analysis 

Beside the cluster lens SDSS1004+4112, whose four bright images are well separated 
from each other (Inada et al. 2003b), the images in the five galaxy-scale lenses are separated 
by only 1-2 arcsec and their PSFs are therefore blended. In each of these systems, we first 
extracted the total counts in each band (full, soft, and hard) within circular apertures 3-4" in 
radius that contain most of the source counts (see Figured]). We estimated the background 
using concentric annular regions with inner and outer radii of ~10" and ~20", respectively. 
We then used PSF fitting to model the relative count rates between the images, fixing the 
relative image positions to the HST positions from the CASTLES website. We binned the X- 
ray images in 0"0984 pixels, and fit these binned images using Sherpa to minimize the Cash 
C-statistic between the observed and model images. The count rates were then normalized 
by partitioning the background-subtracted total source counts based on the relative count 
rates for each image obtained from the PSF fits. Although some images in these five lenses 
are well-separated and their count rates could also be estimated using aperture photometry, 
we find that the count rates from the two methods are consistent and report only the results 
based on the PSF fits. 

In SDSS 1004+4112, the four images are well separated (from ~4" to ~15"), so we 
computed the photon count rates using aperture photometry. Since the foreground lens is a 
cluster, there is non-negligible X-ray emission from the cluster at the positions of the lensed 
images (Ota et al. 2006). We chose a circular region of ~1"5 radius around each image as 
the source region. We estimated the local background in partial annuli at the same radius 
from the cluster center as the quasar images in order to estimate, and then subtract the 
contamination from the lensing cluster. We then corrected for the finite aperture size to 
determine the total count rates in Table HI 

In all six lenses, the count rate of each image was corrected for Galactic absorption 
and absorption in the lens (see Section @J. The final results are presented in Tables [THS1 
We note that the total count rates can differ slightly from the sum of the soft and hard 
band count rates because the PSF fits were done separately for the full, soft, and hard band 
images, and because of the absorption corrections. Finally, we also constructed a stacked 
image (combining all epochs) using the best-fitting image positions obtained from the fits to 
register the epochs and then smoothing with a suitably sized Gaussian kernel. These images 
are shown in Figure [TJ 
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4. Spectral Analysis 

We extracted three sets of spectra from each lens: the spectra of the individual images for 
the individual epochs, the stacked spectra of the individual images combining all the epochs, 
and the stacked spectra of all images. First, we extracted spectra of the individual images us- 
ing the CIAO tool "psextract" in circles of radii about 0"5-0"8 (for SDSS 1004+4112 we used 
1'.'5) centered on the positions from the PSF fits to each observation. Because of the superb 
angular resolution of Chandra, the normal X-ray background contributes little to the source 
spectrum, and the "background" is dominated by contaminating emission from adjacent im- 
ages. To estimate the contamination of image A by image B, we measured the spectrum in 
a circular aperture of the same size at the position X that A would have if reflected through 
image B (i.e., X-B-A with XB = BA). For the cluster lens SDSS 1004+4112, we extracted 
the source spectra using 1'.'5 radius apertures and the background spectra from concentric 
partial rings at the same radius from the cluster center to model the cluster contamination. 
Images C and D of SDSS 0924+0219 are faint and poorly resolved (see Figure [TJ, so we 
extracted a single spectrum from a polygonal region containing images A, C and D. After 
extracting the spectra of the individual images for each epoch, we also obtained combined 
spectra for each image by stacking the spectra from the different epochs and the correspond- 
ing rmf and arf files weighted by their exposure times to increase the S/N ratio of the spectra. 
Finally, we obtained a spectrum combining all the images and all the epochs for each lens. 

We analyzed the spectra using XSPEC V12. We first modeled the spectra of the indi- 
vidual images extracted from the individual epochs using a simple power-law model modified 
by Galactic absorption (Dickey & Lockman 1990) and neutral absorption at the lens over 
the 0.4-8.0 keV band. In the fitting process, we assumed the same fixed Galactic absorption 
for all the images of a lens, but allowed the absorption from the lens galaxy to vary indepen- 
dently. We further assumed that the power-law index Y of the source spectrum was the same 
for all images. In general, this may not be true either because of the intrinsic spectral vari- 
ability modulated by the time delays between the images, or because of energy- dependent 
microlensing, but it is a good average approximation. The spectral fitting results are pre- 
sented in Tables [7l4T2| and Figure [2] shows the evolution of the power-law index V in each 
lens. Besides the power-law index V and the column density Nh of the neutral absorption in 
the lens, we also give the estimated absorption-corrected energy flux for the soft and hard 
bands. 

To better constrain the lens absorption and power-law indices of the quasar, we fit the 
stacked spectra combining the individual epochs for each image. These spectra are shown in 
Figures [3HH1 When fitting these stacked spectra, we also added Gaussian emission lines to 
the model. We detected differential absorption between some of the images in all the lenses 
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except for SDSS 0924+0219 and HE 1104-1805. We detected FeKa emission lines in all 
lenses and the Ni XXVII Ka line in Q 2237+0305. We list the spectral fitting results for 
the stacked spectra in Tables I7HT21 for the photon indices and lens absorption and in Table 
[TBI for the emission lines. After obtaining the best fits to the stacked spectra, we calculated 
the ratios of the count rates between the absorbed and unabsorbed models in the full, soft, 
and hard bands, respectively. These mean ratios were then applied to the absorbed count 
rates of the individual epochs and images to estimate the absorption corrected count rates. 
For SDSS 0924+0219 we use the same absorption correction for images A, C and D. We 
also modeled the stacked spectra allowing different power-law indices between the images. 
Finally, we fit the stacked spectrum combining all images and all epochs for each lens, and 
these results are also presented in Tables I7T-TT21 



4.1. Photon Indices 

We show the time evolution of the power-law index V of each lens in Figure [2J To test 
for spectral variability, we can compare the fits to the individual epochs to their mean and 
its uncertainty (dashed lines in Figure [2]). The results for QJ 0158—4325, HE 0435—1223, 
HE 1104-1805, SDSS 0924+0219, and SDSS 1004+4112 are all consistent with a constant 
spectral index. Only Q 2237+0305 shows clear evidence for spectral variability, with a con- 
stant T model having \ 2 — 38.4 (dof=17), so we detect spectral variability in Q 2237+0305 
at 99.8% confidence. This spectral variability is mainly caused by epochs 15 and 16, where 
the power-law index drops suddenly from ~2.0 to ~1.4. These two epochs contribute ~18.0 
out of the total \ 2 — 38.4, and the spectral variability would not be statistically significant 
without these two points. 

Spectral variability can be either intrinsic variability or caused by chromatic microlens- 
ing between the soft and hard X-ray bands. For example, the spectrum of a microlensed 
quasar image more strongly magnified in the hard X-ray band will appear harder than the 
other images. Distinguishing these two possibilities is particularly important for Q 2237+0305, 
since in this lens we observed more microlensing variability amplitude in the hard X-ray band 
than in the soft. As can be seen from Figures 2 and 3 of Chen et al. (2011), around date 
5400 (epochs 15 and 16), there is no significant microlensing variability in images A, B, and 
D based on the relatively flat B/A, D/A, and B/D flux ratio light curves. Although there 
appears to be some microlensing activity in image C around this time, image C is the faintest 
of the four images at this time and contributes less than 10% of the total flux (Table [121 
fit No. 15, 16), so it cannot be responsible for the sudden decrease in the power-law index. 
Moreover, if we fit the spectra of the four images independently at these two epochs, we find 
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similar photon indices for all images. We conclude that the spectral variability at these two 
epochs is intrinsic, instead of due to chromatic microlensing. We also compare the average 
power-law indices of the individual epochs (dashed lines in Figure with those obtained 
from fitting the stacked spectra combining all the epochs (red points in Figure [2]), and find 
that they are consistent. 

4.2. Absorption in the Lens 

We report our estimates for the amount of absorption in the lens from analyzing the 
combined spectra of each image in Tables ITHT21 (the 5 th column). For absorption we focus 
on the stacked spectra since the changes in the intrinsic spectra are small or non-existent 
and the absorption does not change between epochs. The stacked spectra of QJ 0158—4325, 
SDSS 0924+0219, HE 1104-1805, SDSS 1004+4112, HE 0435-1223, and Q 2237+0305 are 
shown in Figures [3HS1 along with their best-fit models (a power-law modified by Galactic 
and lens absorption plus Gaussian emission lines). Figures I9TTTT1 show the two-parameter 
confidence contours for the estimated column densities Nh and the spectral index T of the 
source. For the two-image lens QJ 0158-4325, we detect differential absorption between 
image A and B (Figure |9]). We detect no lens absorption in HE 1104—1805 (Figure [TOl . For 
SDSS 0924+0219, images C and D are not well separated from image A, so we extracted only 
one spectrum combining images A, C and D, labeled by ACD in Figure [TTJ (bottom right 
panel). We detected no neutral absorption in SDSS 0924+0219. We detected differential 
absorption between images in SDSS 1004+4112, HE 0435-1223, and Q 2237+0305. Dai et 
al. (2003) detected neutral absorption in images A and C of Q 2237+0305, but not in images 
B and D (see also Agol et al. 2009). The latest Chandra data (stacking over 20 epochs with 
a total exposure time over 290 ks) shows that there is neutral absorption in all 4 images, 
where image D suffers the most significant absorption from the foreground lens. 

4.3. Metal Emission Lines 

We detect Fe Ka emission lines in all six lenses. Table [13] reports the rest frame mean 
energy, E line , Gaussian line width, <Ti ine , the equivalent width, EW, and significance for each 
lens. Figure [121 compares the EWs of the individual images, and Figure [TBI shows the average 
(over all images) of the line centers and equivalent widths. We do not detect the Fe Ka line 
for image B of SDSS 0924+0219, QJ 0158-4325, or HE 1104-1805. In HE 0435-1223, we 
detect the FeKa line in images A, B, and C, but not D. In Q 2237+0305, we detect the 
FeKa; line in all four images, confirming the earlier detection of the line in the brightest 
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image A (Dai et al. 2003). In addition, we detect at low significance the Ni XXVII Ka line 
(at E ~ 7.8 keV) in images A and D of Q 2237+0305 (Figure E and Table H3J for the first 
time in a lensed quasar, and we list them in Table [13] for completeness. While adding one 
FeKa line component is sufficient for most spectra, in SDSS 0924+0219, we find that two 
Fe Ka line components can significantly improve the fit for both the spectrum ACD and the 
total spectrum. We detect one blueshifted (5.9 keV) and one redshifted (7.1 keV) FeKa 
line components in both of these spectra with an average line energy of 6.5 keV and total 
EW=1.2 keV. Since we cannot completely resolve the A, C, and D images of this lens, it 
is possible that the two components are from two images. It is also possible that the two 
components are from a single image, possibly from the brightest image A, caused by a caustic 
dissecting the inner accretion disk. In HE 1104—1805, we formally detect the FeKa line at 
the 99% significance level in the spectra of image A and the total image. However, the lines 
are hard to observe by eye, because they are at the energy where the effective area of the 
mirror changes significantly, and it is difficult to determine the continuum by eye. 

The Fe Ka line has a rest-frame energy range of 6.4-6.7 keV depending on the ionization 
state. Most of the FeKa lines we detect have rest-frame energies consistent with the neutral 
line energy at 6.4 keV. In both HE 0435-1223 and Q 2237+0305, we see shifts in the FeKa 
line energy between the images. We detect a redshifted line, 2.5<r from the rest-frame energy, 
in image C of Q 2237+0305 with E = 6.08+0.13 keV. We detect both a redshifted component 
at E = 5.92 keV, 2.5cr from the rest-frame energy, and a blueshifted component at E — 7.12 
keV, > 5a from the rest frame energy, for image A+C+D of SDSS 0924+0219. We also 
detect blue-shifted or ionized lines in image B of HE 0435-1223 with E = 6.72 ± 0.08 keV 
and image D of Q 2237+0305 with E = 6.62 ± 0.06 keV. In HE 1104-1805 we detect a 
blueshifted iron line for image A at E = 6.84 keV, 3a from the rest-frame energy. A small 
fraction of the lines are resolved as broad lines, while most are unresolved, with lower limits 
consistent with narrow lines. The equivalent widths of the lines are significantly larger than 
those detected in unlensed AGN with similar X-ray luminosities, which we will discuss in 
Section 15.31 

5. Microlensing Analysis 

Our microlensing analysis is based on the image count rates presented in Tables dHHJ 
In principle, both microlensing by stars in the foreground galaxy and intrinsic variability 
of the background quasar can cause variability, and intrinsic variability modulated by time 
delays can mimic microlensing in any given eopch. In Q 2237+0305, the time delays are 
observationally constrained to be less than one day (Dai et al. 2003), and the intraday 
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variability of the background quasar is measured to be only ~1% of the total flux (Chartas 
et al. 2001; Dai et al. 2003). The time delays are longer in the other galaxy lenses, (less 
than 15 days in HE 0435-1223, Kochanek et al. 2006; Courbin et al. 2011; about 15 days 
in SDSS 0924+0219, Inada et al. 2003; < 15 days in QJ 0158-4325, Morgan et al. 2008; 
and ~162 days in HE 1104—1805, Morgan et al. 2008). The effects of source variability 
will be most significant in the cluster lens SDSS 1004+4112, where the time delays between 
some of the images are several years (e.g., image C leads A by ~822 days, image D lags 
A by lat least 1250 days, Fohlmeister et al. 2007, 2008). We focus on the evolution of 
the flux ratios between pairs of images to minimize the effects of intrinsic variability. If 
we were focused on flux ratios at a single epoch, the time delay is the relevant time scale 
over which we need to consider the effects of intrinsic variability. This would be a major 
problem for SDSS 1004+4112 with its long A/B to C or D delays. However, when we focus on 
microlensing through the time variability of flux ratios, we can prove mathematically that the 
relevant time scale for being affected by intrinsic variability is the shorter of the time delays 
and the sampling time scale. Hence for images C and D of SDSS 1004+4112, we need only 
worry about the intrinsic variability on times scales of order 90 days rather than the multi- 
year A/B to C/D time delays. Figures [TlHTHl show the evolution of the full, soft, and hard 
band flux ratios for QJ 0158-4325, HE 1104-1805, SDSS 0924+0219, SDSS 1004+4112, 
and HE 0435—1223, respectively. In these figures, we slightly offset the observation dates 
for the soft and hard bands for clarity. We have already discussed Q 2237+0305 in Chen et 
al. (2011). 



5.1. X-ray Microlensing 

In each lens system, we first fit the microlensing light curves (flux ratios) in each of the 
three bands as constants to test for the existence of microlensing variability. For example, 
in the two-image lens QJ 0158-4325 (Figure [T41) . we obtain x 2 — 9.7 (dof=5) with a null 
hypothesis probability of NP = 0.08 in the full band, \ 2 = 6.4 with an NP = 0.27 in 
the soft band, and \ 2 = 8.92 with an NP = 0.11 in the hard band. Thus we weakly 
detect microlensing in the full X-ray band. For HE 1104—1805 (Figure [i~5]) . we obtain 
X 2 = 81.4 (dof=8) with NP = 2.5 x 10~ 14 in the full band, \ 2 = 33.4 with NP = 0.00005 
in the soft band, and \ 2 = 44.9 with NP = 3.8 x 10 -7 in the hard band. Therefore we 
detect microlensing with high confidence in all three bands for HE 1104—1805. In the four 
image lenses SDSS 0924+0219, SDSS 1004+4112, and HE 0435-1223, where we have to 
test the existence of microlensing signal from each flux ratio combination, we summarize 
the microlensing statistics in Tables H3HS3 In SDSS 0924+0219 (Figure [16] and Table ED, 
among the 6 pairs of flux ratios, the B/A pair (images A and B are relatively brighter) gives 
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the best statistical evidence for microlensing with a full band probability of 99.91%. The 
second most significant pair, C/B, shows microlensing at 94% confidence for the full band. 
We also detect microlensing in SDSS 1004+4112 with high confidence (Figure [T7] and Table 
[Tol) . For example, considering the B/D pair in the full band, we have a x 2 — 144.8 for dof = 4 
with a null hypothesis probability NP = 2.6 x 10~ 30 , in the soft band, we have a x 2 = 72.2 
with an NP = 7.7 x 10" 15 , and in the hard band, a \ 2 = 73.1 with an NP = 5.0 x 10~ 15 . 

The microlensing statistics of HE 0435—1223 are given in Table [TBI We observe from 
Figure [TBI that the flux ratios between images C and B are consistent with a constant, with 
X 2 = 0.7, 0.96 and 0.31 (dof =3) for a model with no variability in the full, soft, and hard 
bands. This suggests that the variability of images B and C is mostly intrinsic during the 
observations reported in this paper. This can also be seen from the similarity between the 
light curve of B/A and that of C/A, or from the similarity between the light curves C/D 
and B/D (Figure [TB]) . From Figure [TBI and Table [TBI we see that there is strong evidence 
for microlensing in image A of HE 0435—1223. For example, when fitting B/A ratio as 
a constant (dof=3), we obtain a f = 40.8 with an NP = 7.3 x 10~ 9 in the full band, a 
X 2 = 17.3 with an NP = 0.0006 in the soft band, and a x 2 = 29.8 and an NP = 1.5 x 10~ 6 
in the hard band. As for image D, when fitting the light curves of B/D or C/D pair with 
constants, we find that the light curve is roughly consistent with constant in the full band; 
however, this is not true when fitting the soft or hard band light curves. For example, from 
the B/D light curve, the full band fit gives x 2 — 1-4 (dof=3) with a one sided probability 
of 0.30, whereas the soft band fit gives a x 2 = 6.0 with an NP = 0.11 and the hard band 
fits gives a x 2 = 6.9 and an NP value 0.07. This suggests that image D of HE 0435—1223 
is also gravitationally microlensed and that the lensing is energy-dependent with changes 
in the soft band opposite from those in the hard band (see next section). We assume that 
the variability of images B and C is intrinsic and use their noise- weighted mean, (BC), as 
an estimate of the intrinsic variability. We show the flux ratio {BC) / A between image A 
and (BC) (the intrinsic variability template), and between image D and (BC), (BC)/D, 
in Figure [TBI When we fit (BC) / A as a constant, we obtain x 2 — 101.4 with an NP value 
of 7.9 x 10~ 22 in the full band, x 2 = 33.4 and an NP = 2.7 x 10~ 7 in the soft band, and 
X 2 = 63.0 and an NP = 1.4 x 10 -13 in the hard band. We therefore detect microlensing 
in image A with high confidence. This interpretation is confirmed by the better sampled 
optical light curves presented in Blackburne et al. (2012). Fitting the light curve (BC)/D 
with a constant, we obtain a x 2 = 1-5 with a one sided probability of 0.32 in the full band, 
a x 2 = 9-1 with an NP = 0.03 in the soft band, and a x 2 = 14.1 with an NP = 0.003 in the 
hard band. Therefore, we detect microlensing in image D of HE 0435-1223. 
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5.2. X-ray Energy-Dependent Microlensing 

Next, we test for differences between the soft and the hard energy bands, although only 
the observations of Q 2237+0305 were designed to have enough counts for this purpose. 
For the two-image lenses QJ 0158-4325 and HE 1104-1805, we can simplify compute the y 2 
between the soft and hard band flux ratios. In QJ 0158-4325, we find x 2 — 6-8 (dof=6) with a 
null hypothesis probability NP = 0.35, which suggests no significant chromatic microlensing. 
Among the six observations of QJ 0158—4325 (Figure IT4")) . epochs No. 1, 3, and 6 show 
differences in flux ratios between the soft and hard bands, whereas the other three epochs 
do not. In HE 1104—1805, we find x 2 = 7.2 (dof=9) which also suggests no significant 
chromatic microlensing. 

For the four-image lenses SDSS 0924+0219, SDSS 1004+4112, and HE 0435-1223, we 
construct a x 2 test f° r energy- dependent microlensing using 

( jhaid hard hard r \2 /^soft e soft . soft r \2 

\ u ij b i H'j u ij) \ a ij b i H'j °ij ) \ 

where <i-° ft and cfe ard are the observed count rates in the soft and hard bands at epoch i in 
image j, sf and s^ ard are the separate source light curves in the soft and hard bands at 
epoch i, fi s ° ft and /i^ ard are the mean magnification of the soft and hard bands in image j, 
and 5ij is the microlensing variability at epoch i in image j assumed to be the same for the 
soft and the hard bands. If there is an energy dependence (i.e., 5 4 s ° ft ^ <5^ ard ), a model with 
5^ = d% oft = 5^ rd should fit poorly. For SDSS 1004+4112, we obtain X 2 = 13.4 (dof=5) with 
null hypothesis probability NP = 0.02. For SDSS 0924+0219, we obtain X 2 = 18.6 (dof=7) 
with NP = 0.009. For HE 0435-1223, we find x 2 = 21.6 (dof=3) with NP = 0.00008. 
Therefore, we detect energy- dependent microlensing at 98%, 99.1% and 99.99% confidence 
for these three four-image lenses. 

Although we have detected Fe Ka lines, which probably originate from a different region 
from the hard X-ray continuum and contaminate the microlensing light curves, the typical 
equivalent widths of the line are only 0.2-0.3 keV (observed frame) and do not contribute 
significantly to the hard band X-ray flux. This means that the differential microlensing 
between the soft and hard X-ray bands cannot be attributed to the presence of the iron 
lines. For example, if we use hard band microlensing light curves excluding the spectral 
segment containing the FeKa emission (2.0-2.5 keV observed frame), we find that the Chen 
et al. (2011) results for Q 2237+0305 are unchanged. 

Unfortunately, this statistical detection of an energy dependence is difficult to visualize. 
For example, in Figures fT9TT2~11 we show the RMS flux ratios for various image combinations 
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in SDSS 0924+0219, SDSS 1004+4112, and HE 0435-1223. Except for Q 2237+0305 (Chen 
et al. 2011) and image A of HE 0435—1223, any energy dependence is not visually obvious. 
Recall, however, that these observations are limited to integration times shorter than those 
needed to measure the hard/soft band light curves well. 

5.3. Microlensing of the FeKa Line 

Microlensing of the FeKa line has previously been detected in MG J0414+0534 (Char- 
tas et al. 2002), Q 2237+0305 (Dai et al. 2003), H 1413+117 (Chartas et al. 2007), and 
SDSS 1004+4112 (Ota et al. 2006). For the six lenses with FeKa detections in the com- 
bined spectra (Table [TBI and Figure IT3|) . we can test for microlensing by searching for changes 
in the line (EW or profile) between the images. Unfortunately we can only use the combined 
spectra rather than examining individual epochs because the S/N of the individual epochs 
is too low. This time averaging will smooth and reduce any microlensing signal. 

In principle, if the continuum is microlensed differently from the lines, the equivalent 
widths (EWs) of the lines of the different images will differ. The effect tends to be one-sided 
even though microlensing can both magnify or demagnify a source. The demagnified regions 
tend to be smooth and wide, and thus demagnify large and small regions almost equally, 
while magnified regions have the caustics that will more strongly magnify compact emission 
regions. Thus there will be a tendency for lensed quasars to have either larger or smaller 
EWs than unlensed quasars depending on whether the line emitting region is smaller or larger 
then the continuum. We fit the EWs of the images as constants in each lens system (except 
SDSS 0924+0219, see Figure U3$. For QJ 0158-4325, HE 1104-1805, SDSS 1004+4112, 
and HE 0435-1223, the EWs of the images are well fit by constants of 0.57 keV, 0.38 keV, 
0.60 keV, and 0.88 keV, respectively, whereas the EW of image C of Q 2237+0305 differs 
from those of the other images. For Q 2237+0305, we obtain a best fit constant EW of 0.50 
keV with x 2 — 4.87 (dof = 3) and an NP value of 0.18, which weakly suggests differential 
microlensing effects between the X-ray continua and the FeKa line. Overall, we conclude 
that the flux ratios of the FeKa lines are comparable to those of the X-ray continua. Of 
course this null result does not mean that the Fe Ka line is not microlensed, rather it implies 
that the size of the FeKa emission has to be comparable to the X-ray continuum emitting 
region. 

As another test, we compared the EW-Luminosity relation of these six lenses (Table 
[T71) with those of 88 nearby Seyfert galaxies observed by Suzaku (Fukazawa et al. 2011). 
The inverse correlation between the EW of the iron line and the X-ray luminosity is referred 
to as the X-ray Baldwin effect, and was first discovered by Iwasawa & Taniguchi (1993) 
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for neutral FeKa lines and X-ray luminosities in the 2-10 keV band. Recently, Fukazawa 
et al. (2011) extended the trend to the 10-50 keV band and for ionized FeKa lines. If 
the iron line emission region is smaller than the continuum and therefore more strongly 
microlensed, lenses should on average have higher EWs than unlensed quasars with similar 
luminosities. Figure [22] shows the results where we have divided the local samples into low 
(< 10 22 cm -2 , red circles) and high (> 10 22 cm -2 , black triangles) absorption systems. Our 
six points are marked by blue diamonds, where we have corrected for typical gravitational 
lensing magnifications of 16 for the four image lenses and 8 for the two image lenses. We 
also include the results for MG J0414+0534 (Chartas et al. 2002) and H 1413+117 (Chartas 
et al. 2007). The data from Fukazawa et al. (2011) clearly shows the X-ray Baldwin effect 
of declining EW with luminosity. The gravitationally lensed quasars have systematically 
higher EWs than typical AGN of similar X-ray luminosity (in the 10-50 keV band). To test 
if our data is statistically consistent with typical AGN, we first fit the data in Fukazawa et 
al. (2011) with a power-law model plus intrinsic scatter in equivalent widths, and obtain 

logEW = (6.81 ± 0.51) - (0.21 ± 0.07) logL ± (1.02 ± 0.11) (2) 

with a reduced x 2 = 0.94. Compared to this relation and its scatter, the 8 lenses have 
a x 2 = 25.2 (dof = 8), which is inconsistent with the distribution of typical AGN even 
without considering that all eight lenses lie on one side of the relation. The difference cannot 
be resolved by errors in the magnification correction applied to the luminosity. These are 
uncertain by at most a factor of 2, while a change by a factor of 30 would be needed to 
explain the differences. This suggests that the Fe Ka line is microlensed more strongly than 
the X-ray continua, which in turn implies that the Fe Ka emission region is smaller than the 
X-ray continuum, confirming earlier conclusions (Dai et al. 2003; Popovic et al. 2003, 2006). 



6. Discussion 

In this paper, we report our Chandra monitoring observations of six gravitationally 
lensed quasars: QJ 0158-4325, HE 1104-1805, HE 0435-1223, SDSS 0924+0219, SDSS 1004+4112, 
and Q 2237+0305. We find that simple power-law models plus Gaussian emission lines modi- 
fied by absorption provide good fits to the spectra. We do not detect the predicted reflection 
component, other than the metal emission lines, due to the low S/N in the hard X-ray 
band. We detect intrinsic spectral variability in two epochs for Q 2237+0305, and differ- 
ential absorption between images in QJ 0158-4325, SDSS 1004+4112, HE 0435-1223, and 
Q 2237+0305, but not in HE 1104-1805, or SDSS 0924+0219. From the stacked spectra of 
each lens, we detect FeKa lines in all lenses and in almost all images. In Q 2237+0305, we 
marginally detect the Ni XXVII Ka line in images A and D. 
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We clearly detect X-ray microlensing in all six lenses, and this can now be analyzed to de- 
termine the size of the X-ray emitting regions. The first of these analyses, for HE 0435—1223, 
is presented in Blackburne et al. (2012). Of these 6 lenses, the integration times were de- 
signed to provide accurate measurements in multiple bands only for Q 2237+0305. For this 
lens, we clearly see energy-dependent microlensing. Despite the lower S/N of the observa- 
tions of the remaining 5 lenses, we still succeed in detecting energy-dependent microlens- 
ing in HE 0435-1223, SDSS 1004+4112, and SDSS 0924+0219, but not QJ 0158-4325 or 
HE 1104-1805. 

In general, we do not see shifts in the FeKa energy or EW between images for these 
lenses, except for Q 2237+0305 and HE 0435—1223, which implies that the sizes of the iron 
line emission regions are comparable to those of the X-ray continua. However, since the 
FeKa line is extracted from the time averaged spectrum of each image, any microlensing 
effect has been smoothed. We note, however, that the lenses seem to have significantly higher 
EWs than a sample of 88 nearby Seyfert galaxies observed by Suzaku (Fukazawa et al. 2011). 
This suggests that the FeKa line is more strongly microlensed than the X-ray continuum, 
and this implies a smaller emission region for the iron lines. This comparison is made at 
fixed luminosity, so it is not simply due to comparing local Seyferts with higher redshift 
more luminous quasars. Since the size of the X-ray continuum is already constrained to be 
small, ~10r 9 (Morgan et al. 2008; Chartas et al. 2009; Dai et al. 2010; Blackburne et al. 
2012), the FeKa emission must come from the very inner region of the disk. The existence 
of chromatic microlensing between the soft and hard X-ray bands and the constraints on the 
size of the line emission region have important implications for testing general relativity in 
the strong field regime and the theoretical modeling of quasar accretion disks. 

Although X-ray emission is one of the defining characteristic of AGN, the origin of this 
emission is still unknown. Unlike X-ray binaries, the AGN accretion disk itself is not hot 
enough to generate X-rays. The X-ray emission is generally attributed to inverse Compton 
scattering of soft UV photons from the accretion disk by hot electrons from a "corona" (see 
review of Reynolds & Nowak 2003). Although the origin of the hot electrons in the corona 
is debated, recent microlensing analyses have constrained the extent of this corona to be of 
order 10 r g (Morgan et al. 2008; Chartas et al. 2009; Dai et al. 2010; Blackburne et al. 2012). 
Combining the results of this paper and Chen et al. (2011), we further find that the corona 
might have a temperature gradient with higher temperature electrons occupying a smaller 
volume. 

Several models have been proposed to heat the electrons. In the classic "sandwich" 
corona model (Haardt & Maraschi 1991, 1993) and its extension, the "patchy" corona model 
(Haardt et al. 1994), the gravitational energy of the accretion disk is assumed to be trans- 
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ferred into the disk and corona with constant fractions. This group of models is in conflict 
with the microlensing results, since they predict similar sizes for the X-ray and optical emis- 
sion regions, whereas the microlensing studies uniformly find smaller X-ray emission sizes 
than optical sizes (Pooley et al. 2006; Morgan et al. 2008; Chartas et al. 2009; Dai et al. 2010; 
Blackburne et al. 2011; Blackburne et al. 2012). However, simply requiring that the energy 
transfer process operates only on scales < 10 r g will solve this problem. A second group of 
models assumes that the X-ray emission is composed of a large number of flares triggered by 
magnetic field re- connect ions (e.g., Collin et al. 2003; Czerny et al. 2004; Torricelli-Ciamponi 
et al. 2005; Trzesniewski et al. 2011). The flare models can be considered as realistic re- 
alizations of the "patchy" corona model, and have advantages for interpreting the X-ray 
variability of AGN. Another difference from the steady corona models is that the electron 
distribution is mostly non-thermal in the flare models (Torricelli-Ciamponi et al. 2005). In 
general, the flare models predict similar X-ray and optical emission sizes for the microlens- 
ing analysis, because the corona is composed of a larger number of flares and the flares are 
randomly distributed over the disk. Some recent developments of the flare model specify 
that an avalanche of flares will be triggered by a single flare with a preferred direction of the 
triggered flares to the black hole (Trzesniewski et al. 2011), which will result in a smaller 
X-ray emission size at the disk center. A final set of models heat electrons with shocks. For 
example, Ghisellini et al. (2004) proposed the aborted jet model, where a jet launched near 
the black hole cannot escape the potential well, returns, and collides with other material to 
form shocks. These shocks then convert the kinetic energy of the aborted jet into thermal 
energy. This model predicts compact X-ray emission regions, possibly consistent with the 
microlensing results. 

Given a disk-corona geometry, there should also be a reflection component created by 
the X-ray irradiation of the disk. Complicated radiative transfer processes occur during 
the reflection and the resulting continuum emission peaks at ~ 50 keV and produces metal 
emission lines, such as FeKa (e.g., George & Fabian 1991). Based on the geometry of the 
model, one would expect the reflection region to be much larger than the corona, because 
the reflection can occur anywhere in the disk, a disk wind, or a surrounding torus. With 
the discovery of broad, skewed FeKa lines in many Seyferts (e.g., MCG6-30-15, Tanaka et 
al. 1995), the reflection component near the black hole must contribute significantly to the 
total emissivity. In addition, the reflection components in some Seyferts are measured to be 
constant while the X-ray continuum is variable. Fabian & Vaughan (2003) proposed a light 
bending model, where the majority of the X-ray emission will illuminate the accretion disk 
due to the light bending effect by the black hole, producing the constant reflection compo- 
nent, and the escaping X-ray photons produce the X-ray continuum. If the X-ray source 
moves vertically above the accretion disk, it will produce a variable X-ray continuum and 
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almost constant reflection component. In a competing model, Miller et al. (2008) proposed 
that the reflection occurs in the disk wind and that the redshifted Fe Ka lines are artifacts 
created by a forest of warm absorbers in the X-ray bands. Our microlensing results suggest 
that the Fe Ka emission region is compact and comparable in size or smaller than the X-ray 
continuum region, which prefers the light bending model over the disk-wind model. 

We have composed, for the first time, a modest sample of quasar microlensing light 
curves in the full, soft, and hard X-ray bands, and a sample of Fe Ka measurements in lensed 
quasars. Microlensing of the X-ray emission of lensed quasar is now clearly established in the 
full band (Morgan et al. 2001; Dai et al. 2003; Blackburne et al. 2006; Morgan et al. 2008; 
Chartas et al. 2009; Dai et al. 2010; Chen et al. 2011; Blackburne et al. 2011, 2012), and in 
the FeKa lines (Chartas et al. 2002, 2004; Dai et al. 2003; Popovic et al. 2003, 2006). There 
are well-defined methods of analyzing the results to estimate the sizes of the X-ray emission 
regions (Kochanek 2004; Kochanek et al. 2007; Dai et al. 2010; Blackburne et al. 2012). 
To date, these analyses have yielded only upper bounds on the sizes of the X-ray emission 
regions because of the sparse X-ray light curves and most observations are too short to study 
multiple energy bands or the Fe Ka emission in detail. There are, however, no observational 
or theoretical limits to obtaining the necessary data as long as the Chandra observatory 
continues to operate. Unfortunately, this science will become impossible with Chandra's 
demise because no future mission is designed to have the necessary angular resolution. 
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Table 1. Absorption-Corrected Count Rates for QJ 0158—4325. 



Obsld a 


11556 


11557 


11558 


11559 


11560 


11561 


Date 


2009-11-06 


2010-01-12 


2010-03-10 


2010-05-23 


2010-07-28 


2010-10-06 


Exposure (s) 


5026 


5016 


5035 


4942 


4946 


4944 


Afuii 


28.4±2.9 


33.6±3.3 


34.4±3.0 


35.9±3.3 


37.3±3.5 


30.7±2.9 


A S oft 


13.3±2.1 


18.3±2.4 


20.5±3.5 


20.1±2.5 


18.6±2.5 


20.3±2.4 


Ahard 


13.5±1.8 


14.7±2.2 


13.5±1.8 


15.8±2.0 


17.4±2.0 


10.4±1.5 


Bfuii 


10.7±1.6 


15.2±2.1 


13.3±1.8 


9.7±1.6 


9.0±1.5 


12.2±1.8 


B so ft 


7.7±1.5 


8.5±1.6 


6.2±1.8 


5.0±1.2 


5.4±1.2 


6.5±1.3 


Bhard 


3.9±1.0 


6.5±1.3 


6.7±1.3 


3.9±1.0 


4.1±1.0 


5.4±1.1 


bkgf u ii 


2.97 


1.92 


3.13 


4.97 


3.37 


2.66 


bkg soft 


0.35 


0.35 


0.35 


0.89 


0.18 


0.35 


bkg hard 


2.62 


1.57 


2.79 


4.08 


3.19 


2.30 



a All events have standard ASCA grades of 0,2,3,4,6 and have energies between 0.4-8 keV. We split 
the full X-ray band (0.4-8.0 keV) into a soft X-ray band (0.4-1.3 keV), and a hard X-ray band (1.3-8 
keV). The background counts are extracted from an annulus with inner and outer radii of 5" and 
20", respectively. The units of the image count rates arc 10~ 3 s _1 , and the units of the background 
fluxes are 10 -6 s _1 arcscc -2 . 



Table 2. Absorption-Corrected Count Rates for HE 1104-1805. 



Obsld a 


375 


6918 


6917 


6919 


6920 


6921 


11553 


11554 


11555 


Date 


2000-06-10 


2006-02-16 


2006-03-15 


2006-04-09 


2006-10-31 


2006-11-08 


2010-02-09 


2010-07-12 


2010-12-20 


Exposure (s) 


47425 


4959 


4549 


4865 


5007 


4915 


12756 


12758 


12757 


Af u ii 


25.3±0.8 


8.4±1.5 


10.3±1.7 


7.5±1.4 


16.2±2.2 


13.4±2.3 


10.8±1.1 


23.5±2.1 


11.7±1.1 




15.1±0.8 


3.9±1.0 


5.3±1.3 


3.8±1.1 


9.0±1.9 


7.3±1.6 


5.6±0.9 


11.2±1.2 


5.3±1.0 


Ahard 


10.3±0.7 


4.2±0.9 


5.4±1.2 


4.0±1.1 


5.9±1.4 


6.4±1.3 


4.2±0.6 


11.5±1.4 


5.9±0.8 


Bfull 


13.6±0.6 


21.8±2.8 


17.3±2.5 


16.5±2.3 


8.7±1.5 


13.0±2.3 


14.5±1.3 


9.3±1.1 


14.3±1.2 


Bsoft 


8.6±0.6 


9.3±1.8 


10.4±2.0 


8.9±1.8 


6.1±1.5 


7.3±1.7 


7.7±1.1 


4.6±0.7 


6.4±1.1 


Bhard 


4.5±0.4 


11.3±1.7 


5.8±1.4 


6.5±1.4 


3.6±1.0 


4.8±1.1 


7.0±0.9 


4.2±0.7 


7.2±1.0 


bkgfull 


3.14 


1.99 


2.65 


8.11 


3.94 


2.68 


2.06 


3.18 


3.01 


bkgsoft 


1.02 


0.22 


0.72 


1.80 


0.88 


0.45 


0.09 


0.43 


0.94 


bkghard 


2.13 


1.77 


1.93 


6.31 


3.06 


2.23 


1.98 


2.75 


2.06 



a All events have standard ASCA grades of 0,2,3,4,6 and have energies between 0.4-8 keV. We split the full X-ray band 
(0.4-8.0 keV) into a soft X-ray band (0.4—1.2 keV), and a hard X-ray band (1.2-8 keV). The background counts are 
extracted from an annulus with inner and outer radii of 10" and 20", respectively. The units of the image count rates 
are 10~ 3 s _1 , and the units of the background fluxes are 10 -6 s -1 arcsec -2 . 
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Table 3. Absorption-Corrected Count Rates for SDSS 0924+0219. 



ObsId a 


5604 


11562 


11563 


11564 


11565 


11566 


Date 


2005-02-24 


2010-01-13 


2010-03-10 


2010-05-10 


2010-06-22 


2010-10-05 


Exposure (s) 


17944 


21542 


21264 


21642 


21546 


21543 


Afuii 


6.6±0.8 


12.0±1.3 


8.0±0.8 


18.5±1.2 


13.6±1.2 


9.9±0.8 




3.6±0.6 


7.0±1.0 


3.9±0.6 


9.8±0.9 


7.2±0.9 


5.8±0.6 


Ahard 


2.6±0.4 


4.8±0.6 


3.7±0.5 


7.8±0.7 


5.9±0.9 


3.5±0.5 


Bfuii 


1.7±0.3 


3.2±0.5 


0.9±0.2 


2.8±0.4 


4.6±0.5 


1.9±0.3 


B soft 


1.3±0.3 


1.7±0.4 


0.4±0.1 


1.6±0.3 


2.2±0.4 


1.4±0.3 


Bhard 


0.4±0.1 


1.2±0.3 


0.4±0.1 


1.1±0.2 


2.1±0.5 


0.6±0.2 


Cfull 


0.8±0.3 


0.8±0.3 


1.1±0.3 


1.2±0.3 


0.9±0.3 


1.0±0.3 


Cgoft 


0.5±0.2 


0.8±0.3 


0.7±0.2 


0.5±0.2 


0.5±0.2 


0.7±0.2 


Chard 


0.2±0.1 


0.1±0.1 


0.4±0.2 


0.6±0.2 


0.4±0.2 


0.3±0.1 


Dfull 


0.7±0.3 


2.1±0.5 


0.8±0.3 


1.5±0.3 


1.3±0.4 


1.3±0.3 


D so ft 


0.8±0.3 


1.4±0.4 


0.6±0.3 


1.0±0.3 


1.0±0.3 


0.5±0.2 


Dhard 


0.1±0.1 


0.7±0.2 


0.2±0.1 


0.5±0.2 


0.3±0.3 


1.0±0.3 


bkgf u ii 


1.41 


2.53 


2.57 


2.70 


2.98 


2.35 


bkg soft 


0.22 


0.27 


0.27 


0.72 


0.81 


0.63 


bkg hard 


1.19 


2.26 


2.29 


1.98 


2.17 


1.72 



a We split the full X-ray band (0.4-8.0 keV) into a soft band (0.4-1.2 keV) and a hard band (1.2-8.0 
ke V) . The background counts are extracted from a concentric circle of inner and outer radii 8" and 
15". The units of the image count rates are 10 -3 s _1 , and the units of the background fluxes are 
10 -6 s -1 arcscc- 2 . 
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Table 4. Absorption-Corrected Count Rates for SDSS 1004+4112. 



ObsId a 


5794 


11546 


11547 


11548 


11549 


Date 


2005-01-01 


201 0-03-08 


201 0-06-1 9 


201 0-09-23 


201 1-01-30 

Z. V 1 I V I'M/ 


Exposure (s) 


79987 


5957 


5960 


5963 


5959 


Af u ii 


17.4±0.5 


12.2±1.5 


9.1±1.3 


8.5±1.3 


5.4±1.0 


A so ft 


8.1±0.3 


4.8±1.0 


4.3±0.9 


3.2±0.8 


2.3±0.7 


Ahard 


9.2±0.4 


7.2±1.2 


4.7±0.9 


5.2±1.0 


3.0±0.8 


Bfuii 


22.7±0.6 


16.6±1.8 


9.4±1.3 


12.7±1.5 


8.4±1.3 


B so ft 


11.1±0.4 


8.5±1.3 


3.5±0.8 


4.8±1.0 


2.7±0.7 


Bhard 


11.4±0.4 


7.9±1.2 


5.7±1.0 


7.7±1.2 


5.5±1.0 


Cfull 


18.9±0.5 


12.8±1.6 


18.4±1.9 


12.1±1.5 


21.5±2.0 


Csoft 


9.0±0.4 


5.1±1.0 


8.6±1.3 


5.1±1.0 


9.3±1.4 


Chard 


9.6±0.4 


7.4±1.2 


9.5±1.3 


6.8±1.1 


11.8±1.5 


Dfuii 


11.5±0.4 


17.9±1.9 


12.4±1.6 


9.6±1.4 


16.7±1.9 


D so ft 


5.4±0.3 


10. Oil. 5 


5.9±1.2 


3.9±0.9 


5.7±1.2 


Dhard 


5.8±0.3 


7.7±1.2 


6.2±1.1 


5.3±1.0 


10.1±1.4 



a All events have standard ASCA grades of 0,2,3,4,6 and have energies between 0.4- 
8.0 kcV. We split the full X-ray band (0.4-8.0 kcV) into a soft X-ray band (0.4-1.1 
keV) and a hard X-ray band (1.1-8.0 keV). The units of the image count rates are 
10~ 3 s" 1 . 



Table 5. Absorption-Corrected Count Rates for HE 0435-1223. 



ObsId a 


7761 


11550 


11551 


11552 


Date 


2006-12-17 


2009-12-07 


2010-07-04 


2010-10-29 


Exposure (s) 


9994 


12876 


12757 


12758 


Af u ii 


41.5±2.7 


12.9±1.3 


12.1±1.4 


7.6±1.0 


Asoft 


25.8±2.3 


7.3±1.0 


7.3±1.1 


4.5±0.8 


Ahard 


13.3±1.4 


4.7±0.7 


3.6±0.6 


2.9±0.5 


Bfuii 


10.8±1.2 


12.6±1.2 


8.3±2.0 


6.3±0.8 


Bsoft 


7.5±1.1 


7.7±1.0 


4.3±0.7 


1.9±0.5 


Bhard 


3.4±0.6 


4.6±0.7 


4.0±0.7 


3.6±0.6 


Cfull 


11.5±1.3 


12.3±1.2 


10.5±1.2 


6.7±0.9 


Csoft 


7.5±1.1 


6.4±0.9 


4.0±0.7 


2.2±0.5 


Chard 


4.1±0.7 


5.7±0.7 


5.5±0.9 


4.1±0.6 


Df u ii 


11.5±1.3 


ll.Oil.l 


9.6±2.1 


6.8±0.9 


Dsoft 


5.9±1.0 


7.7±1.0 


5.6±0.8 


3.5±0.7 


Dhard 


5.2±0.8 


3.3±0.5 


4.3±0.7 


2.7±0.5 


bkgfull 


1.48 


2.43 


3.03 


2.38 


bkgsoft 


0.41 


0.83 


0.71 


0.64 


bkghard 


1.07 


1.60 


2.32 


1.74 



a We split the full X-ray band (0.4-8.0 keV) into a soft X-ray band 
(0.4-1.3 keV), and a hard X-ray band (1.3-8.0 kcV). The background 
counts are extracted from an annulus with inner and outer radii of 15" 
and 25", respectively. The units of the image count rates are 10~ 3 s _1 , 
and the units of the background fluxes are 10~ 6 s _1 arcsec~ 2 . 



Table 6. Absorption-Corrected Count Rates for Q 2237+0305. 
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3.6±0.6 


2.9±0.4 


2.8±0.4 


1.1±0.3 


1.2±0.2 


15.3±1.2 


7.0±1.0 


6.0±0.5 


2 


6 


0. 


7 


2 


.0 


XVI 


11537 


55416 


29355 


14.7±1.0 


6.5±0.8 


5.8±0.5 


4.4±0.5 


1.9±0.4 


2.0±0.3 


2.9±0.4 


0.7±0.2 


1.5±0.3 


9.4±0.9 


3.8±0.7 


4.5±0.5 


3 


.3 


0. 


9 


2 


.4 


XVII 


11538 


55471 


29354 


34.0±1.5 


15.3±1.2 


14.7±0.8 


10.3±0.8 


5.2±0.7 


4.1±0.4 


3.3±0.5 


1.8±0.4 


1.2±0.2 


34.3±1.7 


20.7±1.6 


11.7±0.7 


2. 


.7 


0. 


8 


1 


.9 


XVIII 11539 


55524 


9828 


17.6±2.0 


6.0±1.2 


9.6±1.1 


6.4±1.2 


1.1±0.5 


3.3±0.7 


2.8±0.8 


1.0±0.5 


1.1±0.4 


10.0±1.6 


5.8±1.4 


2.9±0.6 


2 


.9 


0. 


7 


2 


.2 


XIX 


13195 


55526 


9830 


16.3±1.9 


6.7±1.4 


7.5±1.0 


5.3±1.0 


1.6±0.7 


2.6±0.6 


2.1±0.7 


1.0±0.5 


0.8±0.4 


8.3±1.4 


3.4±1.2 


3.4±0.7 


3 


.5 


0. 


8 


2 


.7 


XX 


13191 


55528 


9828 


15.7±1.8 


5.7±1.3 


7.6±1.0 


4.4±1.0 


1.8±0.7 


1.9±0.5 


2.0±0.6 


1.1±0.5 


0.7±0.3 


7.2±1.3 


3.0±1.1 


2.9±0.6 


2 


.2 


0. 


7 


1 


.6 



a We split the full X-ray band (0.4-8.0 keV) into a soft X-ray band (0.4-1.3 keV), and a hard X-ray band (1.3-8.0 keV). The background counts are extracted 
from an annulus with inner and outer radii of 10" and 20", respectively. The units of the image count rates are 10~ 3 s -1 , and the units of the background fluxes 
are 10~ 6 s _1 arcscc — 2 . 
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Table 7. Fits to the Chandra spectra of QJ 0158-4325. 



Fit a 


Epoch 


Image 


r 


N H 
xlO 22 

cm -2 


Flux b 

xicr 13 

erg s _1 cm -2 


Flux c 

xl0~ 13 
erg s _1 cm -2 


x 2 » 


P(x» d 


1 
1 


I 
I 


A 
B 


1.86t°;S 


10+ ' 19 

00+ 05 
u - uu -o.oo 


o 7 +0.04 

16+ 001 

U - ±D -0.06 


1 02 +0 ' 20 
1-U -0.22 

33+° 10 

u -°°-0.1() 


1.24(16) 


0.23 


2 
2 


II 
II 


A 
B 


1.95t^ 


00 +U U4 

n 04+ ' 08 

u - u4 -0.04 


0.48_ 07 
n 97 +0.03 


n H7 +U.2U 
54+ ' 13 


1.27(20) 


0.19 


3 
3 


III 
III 


A 
B 


2.001^ 


02+ ' 13 
u - uz -0.02 


0.58±^ 
0.20lH 3 5 


1 05+ uly 
1.UO_q 21 

1fi+ - 10 
u - ao -0.09 


0.82(19) 


0.68 


4 
4 


IV 
IV 


A 
B 


1 7Q+ U17 


00 +U U6 

00+ 09 
u - uu -o.oo 


51+ U - U7 
17 +0.03 


, i 7 +0.23 
± - i '-0.19 

°- 38 -0'.O9 


1.30(19) 


0.17 


5 
5 


V 
V 


A 
B 


L - ,o -0.16 


00 +U U5 
u - uc >-0.03 


51+ U - U7 
u - ol -0.07 

ul ' -0.04 


1 24+ U ' 25 

o 45+ - 12 


1.09(19) 


0.35 


6 
6 


VI 
VI 


A 
B 


2 06 +U - ilJ 


02 +U ' iu 
u - ul -o.oi 


56 +U U3 
21 + 03 


u - yz -0.24 
,0+0.13 
u -°°-0.1() 


1.28(18) 


0.19 


7 
7 


stacked 
stacked 


A 
B 


1 93+ u(J6 


00 +u ui 
n ns+ 03 

u - u,3 -0.02 


53+U.U3 
19+ 001 


99 +U ' UV 
n oq+0.03 


1.11(128) 


0.18 


8 


stacked 


A+B 




00 +(U,i 


85 +(u)4 


i 79 +U.lU 


0.93(129) 


0.70 



a The source and lens redshifts are z s = 1.29 and z; = 0.317, respectively. All models include Galactic 
absorption with a column density of Nh = 0.0195 X 10 22 cm -2 . All derived errors are at 68% confidence. The 
spectral fits were performed within the energy range 0.4-8.0 keV. In the first seven fits, we fit image A and 
B with the same power-law indices but different absorption. In fit No. 8, we extract and stack spectra from 
a circular region of radius ~3" containing both images A and B, and fit it with a simple power-law model 
modified by absorption. 

b Flux is estimated in the 0.4-1.3 keV band (observed 1.3 keV corresponds to 3.0 keV in the rest frame). 
c Flux is estimated in the 1.3-8.0 keV band. 

d P(x 2 /v) is the probability of exceeding x 2 with v degrees of freedom. 
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Table 8. Fits to the Chandra spectra of HE 1104-1805. 



r IT. 


Epoch 


Image 


p 
i 


INlH 

xlO 22 
cm - ^ 


r lux 

xlO" 13 
erg s~~^ cm ^ 


r lux 

xlO" 13 






i 
i 


i 
i 


A 
B 


i s9 +0. 05 
'■■0^-0.05 


n oo+ 001 
u - uu -o.oo 

n nn+ 001 

u - uu -0.00 


0.341°;™ 

o.ist™; 


n 77 +0.05 
n 41+0 03 


0.97(172) 


0.59 


2 
2 


ii 
ii 


A 
B 


1 fiq+ U ' 33 
l.D9_ 26 


°- 28 -CK28 


16+ U ' I)2 

n o 7 +0.03 
U - J ' -0.06 


u -^°_0.18 

1 03+ ' 22 
l.UJ_ 22 


0.60(11) 


0.83 


3 
3 


in 
in 


A 
B 


2 60+ Ud9 


1 23+ iti7 

n 1 ,+ - 25 


51+°- u(i 
45+ - 04 

U.40_Q Q g 


-0.16 
U - ,3Z -0.11 


1.55(7) 


0.15 


4 
4 


IV 
IV 


A 
B 


2 n+ u - 4i 


1 on+ 1,2S 

J-UU_ 1 QQ 

0.00±g;SS 


2Q+ U Ua 

(, OQ+0.00 


049 +0.16 
u - 4z -0 .18 

44+° 17 
u - 44 -0.13 


0.65(6) 


0.69 


5 
5 


V 
V 


A 
B 


9 90+U. 10 
Z - ZO -0.33 


0.28^ 

00+ ' 27 
u - uu -o.oo 


39 +U U3 

21+ 00 
u " 2 -0.08 


44+0.12 
u '^-0.16 

o.24t°; ^ 


1.41(8) 


0.19 


6 
6 


VI 
VI 


A 
B 


1 09+U.34 
1 - SZ -0.24 


U-UU_Q QQ 

n nn +0 19 

u - uu -0.00 


22+°' Ui 
u - zz -0.12 

21+° °° 
u - zi -0.07 


0.49l°;^ 
0.481°;" 


1.13(8) 


0.34 


7 
7 


VII 
VII 


A 
B 


, oo+U.-O 
1 - 88 -0.13 


u - uc> -0.03 

o.oi±°°? 


20+ (, ' l,i 
0-27t°° 2 


42+ U ' 0(i 
u - 4z -0 .07 
-.4-0 (It) 

u - oo -0.08 


0.59(30) 


0.96 


00 00 


VIII 
VIII 


A 
B 


1 7 5+U.14 

'°-0.11 


00 +u ' 14 
u - uu -o.oo 

n nn+ - 07 

u - uu -0.00 


40+ ,J ' ul 
u -* u -o.io 

16+°°° 
U - 1D -0.04 


1 00+"' 12 

41 +°' 07 
u - 4i -0.07 


0.73(38) 


0.89 


9 
9 


vim 
vim 


A 
B 


1 74+U-14 
L - ' 4 -0.13 


00 +Uia 
u - uu -o.oo 

n nn+ ' 05 
u - uu -o.oo 


17 +u.uu 

u - 1 '-0.06 
25 +CU)() 


u '^*-0.08 

64+ - 10 
U - D4 -0.10 


1.10(32) 


0.32 


10 


II 


A+B 


1 QA+U-W 
1 - a *-0.28 


0.38±»f 6 


n 79 +u.u_ 

u -' z -0.06 


1 36 +0 ' 23 
i.ou_q 26 


0.56(14) 


0.90 


11 


III 


A+B 


2 36 +U - 44 
ZOD -0 .23 






69 +U ' 13 

u.oa_ 21 


1.05(10) 


0.40 


12 


IV 


A+B 


ly9 -0 .26 


000 +ui« 


0.45t°;°? 




1.15(9) 


0.32 


13 


V 


A+B 


n o + U.3b 
Z.UO_q 25 


00 +(, ' 2(l 

u - uu -0.00 


50+ (MJ3 


0.82l«;^ 


1.14(12) 


0.32 


14 


VI 


A+B 


1.8b_ 26 


oo +uuy 

u - uu -0.00 


u -*' -0.21 


1 oo +u - 23 

l.UU_ 28 


1.16(11) 


0.31 


15 
15 


stacked 
stacked 


A 
B 


1 - 81 -o.o3 


00 +U U1 
u - uu -o.oo 

n oo+ - 01 

u - uu -0.00 


u - ol -0.03 

n 9 ,+o.oo 
u - zo -o.oi 


u - ' -"--0.03 

52+ 03 
u - oz -0.03 


1.40(48) 


0.03 


16 
16 


stacked 
stacked 


A 
B 


I 76+ u(J& 

II ' D -0.05 
1 Rfi+ - 05 


00 +u ul 
u - uu -o.oo 

n oo +0 - 02 

u - uu -0.00 


30 +(MUI 
n 9 o+0.00 
U ' ZJ -0.01 


74+U W 
n 4q+0-03 


1.38(47) 


0.04 


17 


stacked 


A+B 


1 01+U-U3 
i-»i_0.03 


00 +U U1 

U.UU_Q QQ 


U.O4_ Qg 


1 24 +U.U5 

-■-■^-o.os 


1.01(147) 


0.45 



a The source and lens redshifts are z s = 2.32 and 2; = 0.73, respectively. All models include Galactic 
absorption with a column density of Njj = 0.0462 X 10 22 cm" 2 . All derived errors are at 68% confidence. The 
spectral fits were performed within the energy range 0.4-8.0 keV. In the first nine fits, we fit image A and B 
with same power-law indices but different absorption. In fits no. 10-14 (correspond to epoch no. 2 to 6, in 
which the exposure time is short), we extract the spectrum from a circular region of radius ~4" containing 
both images A and B, and fit it with a simple power-law model modified by absorption. In fit no. 17, we 
extract and stack the spectra from a circular region containing both images A and B, and fit it with same 
model. 

b Flux is estimated in the 0.4-1.2 keV band (observed 1.2 keV corresponds to 4.0 keV in the rest frame). 
c Flux is estimated in the 1.2-8.0 keV band. 

d P(x 2 /f) is the probability of exceeding \ 2 with v degrees of freedom. 
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Table 9. Fits to the Chandra spectra of SDSS 0924+0219. 



Fit a 


Epoch 


Image 


r 


N H 

xlO 22 
cm 


Flux b 

xlO" 13 
erg s _1 cm -2 


Flux c 
xl0~ 13 
erg s _1 cm - 2 




P( X 2 /vf 


1 
1 


I 
I 


A+C+D 
B 


2 <57+ - 52 
z -°'-0.44 


13+ - 20 
u - iJ -0.13 

00+ 09 

u - uu -0.00 


12+ 001 
u - iz -0.03 

0.021°;°° 


n 1 c+0.05 
°- 16 -0.05 

02+ ' 01 
u - uz -o.oi 


1.23(11) 


0.26 


2 
2 


II 
II 


A+C+D 
B 


^• aa -0.13 


u ' uu -0.00 

00+ ' 08 
u - uu -o.oo 


26 +U ' 03 

u - ZD -0.03 

05 +001 

U - UO -0.01 


3l+°- U4 

u -° -0.04 
n nc+0-01 

0.06I 001 


0.67(30) 


0.92 


3 
3 


III 
III 


A+C+D 
B 


9 9 o+U.41i 
Z - ZO -0.37 


0.l3^;f 3 
0-39±S:|g 


u ' 1J -0.02 
01 +0 '°° 

u - ul -0.00 


U - z °-0.06 
04+ ' 02 


1.09(17) 


0.35 


4 
4 


IV 
IV 


A+C+D 
B 


2 02 +U - i3 
z - uz -0.11 


00 +lu)4 
u - uu -o.oo 

00+ 05 

u - uu -0.00 


O.sot^ 

o.o5±g : gi 


62 +(, ' l,V 

o.iotg;g 


0.79(44) 


0.84 


5 
5 


V 
V 


A+C+D 
B 


z - i8 -0.13 


00 +U U4 

u ' uu -0.00 

03+ 11 

U - UJ -0.03 


n 90+0.02 

u - ZJ -0.02 

06+ 01 
u - UD -o.oi 


0.381°°, 
10+ 002 


0.89(34) 


0.65 


6 
6 


VI 
VI 


A+C+D 
B 


2 07+ u i!1 
z - u ' -0.14 


00 +U ' (,S 
u - uu -o.oo 

u - uu -o.oo 


19+ UU2 
u ' iM -0.02 

03+ ' 01 


o.36^-.S 

0.061°;°! 


0.95(27) 


0.53 


7 
7 


stacked 
stacked 


A+C+D 
B 


9 19+U.UV 

z - lz -0.05 


00+ U U3 
u - uu -o.oo 

n 02+ ' 03 

u - uz -0.02 


21 +tU,i 

u - zl -o.oi 

04+ ' 00 
u - U4 -o.oo 


3<5+ (, - l)2 

06+ ' 01 
u - UD -o.oi 


0.90(161) 


0.81 


8 
9 


stacked 
stacked 


A+C+D 
B 


9 90+U.UV 

z - zo -0.05 

2 14+ - 15 
z - i4 -0.14 


00 +tU)4 
u - uu -0.00 

00+ 07 

u - uu -0.00 


21+ tu,,) 
u - z± -0.02 

04+ 00 
u - U4 -o.oi 


0.361°;^ 
0.06±g;gi 


0.91(125) 
0.69(27) 


0.77 
0.89 


10 


stacked 


A+B+C+D 


2 19 +u.uy 

z - la -0.08 


02 +U ' U3 
u - uz -0.02 


97+ u - ul 
UiZI -0.02 


0.49±li;!! 3 


0.95(166) 


0.65 



a The source and lens redshifts are z s = 1.524 and z; = 0.39, respectively. All models include Galactic absorption 
with a column density of Njj = 0.0375 X 10 22 cm -2 . All derived errors are at 68% confidence level. The spectral fits 
were performed within the energy range 0.4—8.0 keV. In the first seven fits, we fit images A+C+D and B with the 
same power-law indices but different absorption. In fits 8 and 9, we fit the stacked spectra of A+C+D, and of B 
separately. In fit No. 10, we extract and stack spectra from a circular region of radius ~3" containing image A, B, 
C and D, and fit it with a simple power-law model modified by absorption. 

b Flux is estimated in the 0.4-1.2 keV band (observed 1.2 keV corresponds to 3.0 keV in the rest frame). 

c Flux is estimated in the 1.2—8.0 keV band. 

d P(X 2 / l/ ) is the probability of exceeding \ 2 with u degrees of freedom. 
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Table 10. Fits to the the Chandra spectra of SDSS 1004+4112. 



Fit a 


Epoch 


Image 


r 


N H 

xlO 22 
cm -2 


Flux b 

xlO" 13 
erg s _1 cm~ 2 


Flux c 
xl0~ 13 
erg s _1 cm -2 




P( X 2 /v) d 


1 
1 
1 
1 


I 

I 
I 
I 


A 
B 
C 
D 




u ' u,;> -0.03 
05+ 01 

10+ 002 

14+ 002 


n 07+0.01 
u ' z '-0.01 
34+ 01 

27+ - 01 
u - z '-o.oi 

14+ - 01 

U- 1 -0.01 


0.82±°;°§ 
70+ 03 

u - ' u -0.03 
40+ ' 02 


1.13(559) 


0.02 


2 
2 
2 
2 


II 
II 
II 
II 


A 
B 
C 
D 


2.08t u ;i; 


n oo+U.bU 
u " 3 " 3 -0.33 

n ns+ 09 

00+ 03 
u - uu -o.oo 

n n7 +0.07 


u ' lo -0.06 
n 97 +0.03 
u - z '-0.04 
24+ ' 01 
n 27 +0.04 


0.5235:3 

0.591°;!? 
n 44+0-09 

u - oa -o.io 


0.95(34) 


0.54 


3 
3 
3 
3 


III 
III 
III 
III 


A 
B 
C 
D 


2 00 +u iK 


00 +u i " 
u - uu -o.oo 

00 +0 ' 07 
u - uu -o.oo 

01+ 06 
u - u± -o.oi 


u - 1D -0.03 
1fi+ 002 

S1+ 04 
n iq+0.03 
°- 18 -0.03 


n oo+U.UO 

u-jj-o.oe 

no + 0.07 

n cc+0.11 

0.66_ 10 
46+ ' 08 


1.05(29) 


0.40 


4 
4 
4 
4 


IV 
IV 
IV 
IV 


A 
B 
C 
D 


L - '*-0.17 


0-59^0.59 

u ' uu -0.00 

00+ 02 
u ' uu -0.00 

14+ - 23 


n n7 +U.U5 
u - u '-0.04 

u ' zu -0.05 
20+ 01 

12+ 002 


5fi+ UiV 
U.OD_ 20 

62 +0 ' 17 
u - oz -0.14 

63+° 13 

47+ - 12 
u -^' -0.11 


0.98(24) 


0.48 


5 
5 
5 
5 


V 
V 
V 
V 


A 
B 
C 
D 


1 Q7+U-14 
1 - a '-0.18 


00+ U UU 
u - uu -o.oo 

00+ 05 

u - uu -o.oo 

00+ 02 

u - uu -o.oo 

1Q+° 16 


10+ tu)2 

11+ 002 
u - 11 -0.02 
o 7 +0.04 
u -' 5 ' -0.04 

21 +0 ' 03 
UiZ± -0.03 


n 9 ,+u.U6 

u - zo -0.05 
67+ ' 11 

u - D '-o.io 


1.02(29) 


0.43 


6 
6 
6 
6 


stacked 
stacked 
stacked 
stacked 


A 
B 
C 
D 


i - JJ -0.03 


o.o4l" o: ^ 

0.09t° ;°i 
0.09±°;81 
1S+ 002 

u - io -0.()2 


24+ u ul 
u ' z4 -0.01 

30+ 01 

27+ ' 01 
u - zl -0.01 

16+ 001 
u - io -0.()l 


u - oo -0.02 
„ --+0.03 
u -' '-0.03 

69+ ' 02 
u - oa -0.02 

4S+ 02 


1.15(659) 


0.01 


7 
8 
9 
10 


stacked 
stacked 
stacked 
stacked 


A 
B 
C 
D 


i- y o_ .05 

1 qs+005 
, qo+0.05 
-, o 7 +0.10 


00+ U U4 
u - uu -o.oo 

00+ 02 

u - uu -o.oo 

01+ 04 
u ' ul -0.01 

06+ 06 


24+ tu,i 
U - Z41 _0.01 

Q OO + 0.01 

29+ 01 
u ' za -0.01 
n ic+0.01 
°- 16 -0.01 


u - oy -0.03 
n 7^+0.03 
u -'°-0.03 
n 71 +0.03 
u - ' -0.03 
n 47+0.03 
u -^'-0.03 


1.11(149) 
0.84(180) 
0.73(167) 
0.89(116) 


0.17 
0.94 
1.00 
0.80 


11 


stacked 


A+B+C+D 


1 - ao — 0.04 


04 +U ' U2 


1 ' i4 -0.02 


2 8Q+ (U,V 
z - oa -0.07 


0.98(375) 


0.61 



a The source and lens redshifts are z s = 1.734 and z\ = 0.68, respectively. All models include Galactic absorption 
with a column density of Njj = 0.0111 X 10 22 cm -2 . All derived errors are at 68% confidence level. The spectral fits 
were performed within the energy range 0.4-8.0 keV. In the first 6 fits, we fit image A, B, C and D simultaneously 
with the same power-law indices but different absorption. In fits 7-10, we fit the stacked images of A, B, C and D 
separately. In fit 11, we extract and stack spectra from a circular region of radius ~10" containing image A, B, C 
and D, and fit it with a simple power-law model modified by absorption. 

b Flux is estimated in the 0.4-1.1 keV band, (observed 1.1 keV corresponds to 3.0 keV in the rest frame). 

c Flux is estimated in the 1.1-8.0 keV band. 

d P(x 2 /^) is the probability of exceeding \ 2 with u degrees of freedom. 
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Table 11. Fits to the Chandra spectra of HE 0435-1223. 



Fit a 


Epoch 


Image 


r 


N H 

-i .-,99 

xlO 22 
cm -2 


Flux b 
xlO -13 
erg s _1 cm~ 2 


Flux c 
xlO 13 

erg s _1 cm -2 


x2M 


P( X » d 


1 
1 
1 
1 


I 

I 
I 
I 


A 
13 
C 
D 


2 00+° 11 
z - uu -o.io 


0-08_ 08 
0.00+°° 3 

u w -0.()(] 

00+ 03 
u ' uu -0.00 

12+° 16 


41+ 04 

0.18+ ;' , 2 

v u — 0.02 

17+ 002 

u - i '-0.02 

14+ 002 


95+ ' 11 
u - 9o _0.10 

0.35 +0 n °l 

u ^-0.05 
34+0 05 

o 7 +0.06 


1.02(58) 


0.43 


2 
2 
2 
2 


II 
11 
11 
II 


A 
13 
C 
D 


1 Q 3 + Uil 


00 +U ' 12 
u - uu -o.oo 

o oo+SSS 

u,vu -0 .[)() 

o.oo+2-SS 

— UAH) 

o oo +004 

u - uu -o.oo 


17+ U - U2 

u -l'-0.02 

18+SS2 

u "-0.02 
n 1^+0.02 
0.16 n no 
— U.Uz 

u - la -0.02 


°- 36 -0.05 



— 0.05 
n o^+0.06 
0.OD n 

— U.Uo 
n 9 o+0.04 


0.96(50) 


0.56 


3 
3 
3 
3 


III 
III 
III 
III 


A 
B 
C 
D 


L - ' 1 -0.11 


o oo +um 

u - uu -0.00 

°- 06 -0.06 

00+ 11 
u - uu -o.oo 

00+ 09 

u - uu -0.00 


12 +U ' U2 
u ' iz -0.02 

09+ 01 
12+ 002 
12+ 002 


°- 36 -0.05 
n , 9 +0.05 
u - oz -0.05 
n q7 +0.06 
U - J ' -0.05 
n o fi +0.06 
U - JD -0.05 


1.15(39) 


0.24 


4 
4 
4 
4 


IV 
IV 
IV 
IV 


A 
B 
C 
D 


1 58+ Ui!7 
l.OO_ 23 


o oo +u ^ 

u - uu -o.oo 

n 31 +0.55 
u -° -0.29 
00+ 08 

u - uu -o.oo 
n ns+01 9 


u ' ua -0.02 

n 05 +0 01 

07+ 01 
u - u '-o.oi 

0.06I 001 


n ,,+u.U6 

u - J,:> -0.06 
34+0-07 
u - o ^-0.06 

24+ ' 06 
U -^-0.05 
n 97 +0.06 
u - z ' -0.05 


1.35(23) 


0.12 


5 
5 
5 
5 


stacked 
stacked 
stacked 
stacked 


A 
B 
C 
D 


1 R1+U.U5 

J-si-o.os 


00+ U U3 
u - uu -o.oo 

n oo+ 003 
u - uu -o.oo 

n nn+ 03 

u - uu -0.00 

0.02t°;°J 


n is+u- ()i 
u - i8 -o.oi 

12+ 001 
u -l-'_0.01 

n 13+001 
u - 11 -o.oi 


47+0-03 
-0.03 
n 09+O.O2 
u - oz -0.02 
33+0-02 

3i + 02 

u - ol -0.02 


1.09(187) 


0.20 


6 
7 

8 
9 


stacked 
stacked 
stacked 
stacked 


A 
B 
C 
D 


2.06tr y 9 

1.871°;" 
1 59+0-11 

1 89+ 011 


0.081H? 
nn+ 05 

u ' uu -0.00 

00+ 04 
u - uu -o.oo 

07+° 12 

u - u '-0.07 


n is+u- ()i 

U - iS -0.02 

13+ 001 

11+ 001 

11+ - 01 
U - L1 -0.01 


42+ (u)4 
n Qo+0.02 
u -° -0.02 
42+ - 05 

30+ ' 04 
U - JU -0.04 


1.04(56) 
1.12(35) 
0.99(41) 
0.94(37) 


0.39 
0.28 
0.49 
0.58 


10 


stacked 


A+B+C+D 


1 90 +UU7 


04 +u (J4 

u - u ^-0.03 


62 +U ' ,)2 

u - DZ -0.03 


1 61+ ,U)8 

1 - D1 -0.08 


0.97(171) 


0.61 



a The source and lens redshifts are z s = 1.689 and z\ = 0.46, respectively. All models include Galactic absorption 
with a column density of Njj = 0.0511 X 10 22 cm -2 . All derived errors are at 68% confidence level. The spectral fits 
were performed within the energy range 0.4-8.0 keV. In the first 5 fits, we fit image A, B, C and D jointly with the 
same power-law indices but different absorption. In fits 6-9, we fit the stacked images of A, B, C and D separately. 
In fit 10, we extract and stack spectra from a circular region of radius ~4" containing image A, B, C and D, and fit 
it with a simple power-law model modified by absorption. 

b Flux is estimated in the 0.4-1.3 keV band, (observed 1.3 keV corresponds to 3.5 keV in the rest frame). 

c Flux is estimated in the 1.3—8.0 keV band. 

d P{x 2 l v ) i s the probability of exceeding \ 2 with v degrees of freedom. 
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Table 12. Fits to the Chandra spectra of Q 2237+0305. 



Fit a 


Epoch 


Image 


r 


N H 

xlO 22 
cm 


Flux b 

xl0~ 13 
erg s _1 cm -2 


Flux c 
xlO" 13 
erg s _1 cm -2 




P(x» d 


1 
1 

1 
1 
1 


T 
1 

I 

I 

I 


A 
1 V 

B 
C 
D 


1 «S+ 007 
1-88-0.07 


n n9 +o.02 
U - Uz -0.02 

00+ 02 
u - uu -o.oo 

o 05+ 002 

u - uo -0.02 
02+ 02 


n t-n+0.02 
u - oz -0.02 

10+ 000 

u - iu -0.02 

20+ 01 
u - zu -o.oi 

11+ 001 

u - il -o.oi 


i o 7 +0.08 
L -°< -0.08 

24+ 03 
u - z4 -0.03 

o.57l°;°i 
o.28t°;° | 


1 .\Jo\ lUo ) 


u.oy 


Z 

2 
2 
2 


TT 
11 

II 
II 
II 


A 

1\ 

B 
C 
D 


! 77+U.lU 
l -> '-0.10 


U - Uz -0.02 

00+ 04 
u - uu -o.oo 

n ns+ ' 08 
02+ 04 

U - UZ -0.02 


"■^"-0.03 

09+ 01 
u ' uy -0.01 

12+ 002 
u - iz -0.02 

12+ 002 

U - iZ -0.02 


1 ns+ u ' 13 

2 t 5+ ' 05 
u - zo -0.05 

45+ - 07 

n oc+0.06 

U - JO -0.06 


l.lO^OU ) 


U.zo 


<? 
O 

3 
3 
3 


TTT 
111 

III 

III 

III 


A 

B 
C 
D 


1 7Q+U-aU 

1 - ' y -0.21 


U - Ub -0.06 

o oo+ - 04 

u - uu _0.00 

o oo +0 - 10 
u - uu -o.oo 

02+° 14 

U - UZ -0.02 


U ' ib -0.03 

n is+ 02 

u - iJ -0.02 
05+ - 01 

07+ 02 


^8+ U - i4 
0.b»_ .12 
34+0.10 

14+0-05 
u - ±4 -0.04 
21+ 06 


1 ao/iy 1 ! 

l.OZ^l 1 J 


n 09 
u.uz 


A 
^± 

4 
4 
4 


TV 
1 V 

IV 
IV 
IV 


A 
B 
C 
D 


1 - 9U -0.15 


n rK+ u - iu 

U.0b_ 05 

00+° 10 
u - uu -o.oo 

00+ 02 
u - uu -o.oo 


21 + ' 01 
u - zl -0.05 

17+ 002 
u - i '-0.03 

U-10_o.05 


1 n7+°- ib 

i - U '-0.16 

50+° 11 

n 47 +0.08 
U-*' -0.09 
n QK+0.08 
°- 36 -0.08 


O Q1 ( A9>\ 


U-DO 


5 
5 
5 
5 


V 

V 
V 
V 


A 

B 
C 
D 


9 n7 +0.39 

z - U '-0.35 


n ifi+ u - i4 

U - ib -0.12 

07+° 16 

u - zl -0.13 
16+ ' 27 


n is+0-03 

09+ 01 
u.uy_ 02 

n n? +0.02 

05+ ' 01 
u - uo -0.02 


D ^q+ u - iu 
U-Oy_o.l3 

21 + 08 
U-^l_o.07 
n 9q +0.06 

u-^y_ .08 

n ic+0.08 

0.16_„ 07 


n ail 1 a\ 

u.y ( 


u.4y 


6 
6 
6 
6 


VT 
V 1 

VI 
VI 
VI 


A 
B 
C 
D 


9 n o+U.i!U 
2-08_ 19 


n n«+ u ' uv 

n ipc+0.10 
u - lo _0.09 
25+° 13 

22+ ' 13 
u - z -o.io 


U-45_ oe 

17+003 

Uii ' -0.03 

09+ 02 
u - uy -0.02 

10+ 002 
u - iu -0.02 


i io+0.14 
1 - li -0.15 
56+ ' 09 
41+0.08 

40+ 08 
u -* u -0.08 


n 07/ x>n\ 
U.y f i^OU j 


U.Oo 


7 
7 
7 
7 


VTT 
V 11 

VII 
VII 
VII 


B 
C 
D 


1.9U_ 20 


u - zy -0.11 
00+ 04 

u - uu -o.oo 
o 1 1 +o io 

u - iJ --0.08 
u - lo -0.09 


u.oo_ 05 

o i^+o oo 

U.lO_ 05 

u.uy_ 02 
n+o ° 2 

u - il -0.02 


z - UJ -0.28 
n oc+0.09 
U.Od_ 09 

Q no + 0.07 

u.oo_ 08 

o 41+00 9 






00 00 00 00 


VIII 
VIII 
VIII 
VIII 


A 
B 
C 
D 


1 64 +U ' 2i 


05 +0 ' U4 

00+ 08 
u - uu -o.oo 

16+° 14 

u ' io -0.11 

27+ ' 18 
u - z '-0.13 


27 +0 ' U2 
u - z ' -0.05 

u.uy_ 04 

06+ 02 
u - UD -0.02 

06+ 02 
u - uo -0.02 


^^^^^-o.is 
si+0 7 

u - r>1 -0.08 
n 07+0. 08 
u -°< -0.09 

0.51±°i| 


0.79(31) 


0.78 


9 
9 
9 
9 


vim 
vim 
vim 
vim 


A 
B 
C 
D 


1- ' ' -0.24 


05 +U ' 12 

o n5+ - 09 

U.UO_ 05 

o io+ - 17 
u - io -o.io 


n os+0.02 
U-^S_ .06 

u - lo -0.03 

0.06t°;» 2 
O.ost™ 2 


0-98±»;ig 
45+0- 11 

25+ ' 08 
u - zo -0.08 

0.39t°;° l 


0.71(29) 


0.88 


10 
10 
10 
10 


X 
X 
X 
X 


A 
B 
C 
D 


1 - 8 '-0.16 


u.uu_ oo 
0.28t°;£ 


0.28tr o4 
0.09t°;g 2 

0.08t°;° 2 
0.05l°;» 2 


U-OU_ 07 

25+ 06 
u - zo -0.06 

29+ ' 09 
u - zy -0.09 


1.27(26) 


0.16 
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Table 12— Continued 



Fit a 


Epoch 


Image 


r 


N H 
xlO 22 
cm -2 


Flux b 
xlO- 13 
erg s -1 cm -2 


Flux c 

xlO- 13 
erg s _1 cm- 2 


x5M 


P( X 2 /v) d 


11 
11 
11 
11 


XI 
XI 
XI 
XI 


A 
B 
C 
D 


2 02+ ' 13 

z ' uz -0.13 


n i4+ - 04 

u ' 14 -0.04 

n m +0 ' 06 

n ns+ - 08 

U - US -0.06 
u - uo -0.04 


n Q9+0 07 
U - yz -0.07 

S2+ ' 04 
U - JZ -0.04 
n 24 +0.03 
u - z4 -0.03 

0.29t°° 04 


3 Ofi+0-25 

J.UD_q 2g 

o.93±g:iJ 

0.641°;}° 

71+0- 10 
u - ' 1 -0.11 


1.08(111) 


0.27 


12 
12 
12 
12 


XII 
XII 
XII 
XII 


A 
B 
C 
D 


2 04 +(, ' i(i 
z - u *-0.09 


03+ ,UI4 

00+ ' 03 
u - uu -o.oo 

00+ ' 08 
u - uu -o.oo 

ii+ - 09 


n 33+0.01 
u - oo -0.07 

3n +0 '00 
u - JU -0.07 

iq+0- 03 
u - la -0.03 


1 85 +l, ' i!) 

62+ 09 
U ' DZ -0.09 
n .-+0.09 
u -°'-0.09 

5fi+ 009 
U - JD -0.09 


0.93(86) 


0.66 


13 
13 
13 
13 


XIII 
XIII 
XIII 
XIII 


A 
B 
C 
D 




10+ ,UI2 
n ,,+0.03 
u - ±,3 -0.03 

n 1S+ ' 07 

U-l»_o.06 

n 14+003 

u - 14 -0.03 


n «q+ u - u4 
U - 89 -0.04 

0.24l°;° 2 

om±° ° \ 

46+ 02 

u -^°-0.03 


09 +U.lU 
z - zz -0.10 
R7+0-05 
u - D '-0.05 

26+ 03 
u - ZD -0.03 

1 34+0.07 


0.86(322) 


0.96 


14 
14 
14 
14 


XIV 
XIV 
XIV 
XIV 


A 
B 
C 
D 


2.00i^ 7 


22 +(J - U7 
u - zz -0.06 

n 16+° 10 

u - MZ -0.38 
u - JU -0.09 


13+ uui 
u - ic> -0.01 

04+ 01 
u - u ^-o.oi 

u - ui -o.oo 

06+ ' 01 

u - UD -o.oi 


U - OD -0.06 

15+0-0 3 

u - io -0.03 

14+ 003 
u ' 1 ^-0.04 

31+0-0 4 

u - ol -0.04 


0.82(66) 


0.85 


15 
15 
15 
15 


XV 
XV 
XV 
XV 


A 
B 
C 
D 


1 40+ ' 18 


12+ ,)42 

u - ±z -0.09 

u - u4 -0.04 
06 +0.21 
u - UD -0.06 

02+ ' 06 
u - uz -0.02 


n+ u - Ui 

u - 4i -o.oi 
n 04+ ' 00 

U - U4 -0.01 

02 +0 '°° 
u - uz -o.oo 

0.07t°;°l 


0-83_ 09 

24+ ' 03 
u ' z4: -0.04 

10+ 002 
u ' lu -0.03 

n 00+0. 05 


1.09(71) 


0.29 


16 
16 
16 
16 


XVI 
XVI 
XVI 
XVI 


A 
B 
C 
D 


1 43+uiy 


06+ ,UIS 

n is+ ' 20 

u - lc> -0.13 
45+0-30 
u -*°-0.22 

u - uo -0.05 


0.08tr o i 

n 02+ ' 00 
u - uz -o.oi 

n oi + 000 
u - ui -o.oo 

04+ ' 00 
u - U4 -o.oi 


n 48+0- 0(i 
0-48_ 06 

15+00 3 

u - ld -0.03 

12+ 003 
u - iz -0.03 

25+ 04 
u - zo -0.04 


0.92(46) 


0.63 


17 
17 
17 
17 


XVII 
XVII 
XVII 
XVII 


A 
B 
C 
D 


2 07+"' 12 
z - u ' -0.12 


U" 3 -1-0. 06 

n 22+ ' 10 

u - zz -0.08 

n nq+ - 12 

u - uy -0.08 
U - zl -0.05 


U - 18 -0.02 

u - uo -o.oi 
0.031°;°° 

0.18t°° 2 


qi+0.0« 
u - yi -0.07 
n oo +0.03 
u ' zz -0.03 

07+ ' 01 
u - u '-0.02 

fiq+ 06 


0.90(120) 


0.78 


18 
19 
20 


XVIII 
XVIIII 

XX 


A+B+C+D 
A+B+C+D 
A+B+C+D 


z - z °-0.30 

1 7O + 0.30 

1 -'°-0.28 
1 ^q+0-39 
l.O»_ 27 


59+ ' 24 
u ' oy -0.18 

n Q1+0.18 
si+ ' 20 

U - Ji -0.17 


0.17+°^ 
0-17tg°oI 

o.i5tg;°oi 


, 97 +U.12 
i - z ' -0.22 
-1 Qo+0. 14 
i " 3z -0.24 
1.44+°i| 


1.37(26) 
1.08(22) 
0.63(20) 


0.10 
0.35 
0.89 


21 
21 
21 
21 


stacked 
stacked 
stacked 
stacked 


A 
B 
C 
D 


±.eo_ 21 


n o C +0.14 

U - OD -0.19 

17+ - 22 
u - ± ' -0.14 

21+ ' 12 
u - zl -o.io 


0.07+°°} 

02+°°° 
u - uz -o.oo 

01+0 00 
u - ul -o.oo 

0-04±g-gi 


u - o,3 -0.07 
n 1Q +0.04 
0.19_ 04 

n7 +0 02 
u - u ' -0.02 

n 1Q +0.03 
U - iy -0.04 


0.72(45) 


0.92 


22 
23 
24 
25 


stacked 
stacked 
stacked 
stacked 


A 
B 
C 
D 


1 qn+ u - u4 
1-9U_ 04 

1 qo+0.08 

2 n4+ ' 10 

z - u4 _0 .09 


0.08tg;g 
n 07+0- 03 

u - u ' -0.03 

n i4+o- 02 

u - i4 -0.02 


°- 64 -0.01 

n lq+o oo 
u - iy -o.oo 
n i4+o-oo 

u - 14 _0 .00 

3n+ 00 

u - JU -0.00 


1 m+ iim 

1 - Li '-0.03 
n oo+O.Ol 

u -°°-o.oi 

23+ 01 

u - zo -o.oi 

52+ 02 
u - oz -0.O2 


1.30(98) 
1.17(97) 
0.99(75) 
0.92(128) 


0.02 
0.12 
0.49 
0.73 


26 
26 
26 
26 


stacked 
stacked 
stacked 
stacked 


A 
B 
C 
D 


1-89-0.03 


07+ ,UU 
u - u '-o.oi 

n nc+0-02 

0.06I .02 

n 03+ 001 

u - 1,5 -0.01 


65 +U ' Ui 
u - DO -o.oi 

U-18_o.oo 
10+O.OO 
u - ±o -0.00 
n oq+0-OO 
u - zy -0.01 


1 iq+0-02 

1 - iy -0.02 
n 00+O.O2 
u - c>z -0.02 
23+ ' 01 

u - zo -o.oi 

52+0' 01 
u - oz -0.02 


1.35(423) 


1.96 x 10" 6 
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Table 12 — Continued 



Fit a 


Epoch 


Image T 


N H 
xlO 22 
cm~ 2 


Flux b Flux c 
xl0~ 13 xlO" 13 
erg s _1 cm~ 2 erg s _1 cm -2 




P( X 2 /v) d 


27 


stacked 


A+B+C+D 1.85±° ° l 


07+ 01 
u - u '-o.oi 


94+O O 1 o 26+ 05 
u - a -0.01 ,3 - ZD -0.05 


1.64(109) 


0.00 



a The source and lens redshifts arc z s = 1.69 and z; = 0.0395, respectively. All models include Galactic absorption 
with a column density of Nh = 0.0551 X 10 22 cm" 2 . All derived errors are at 68% confidence level. The spectral fits 
were performed within the energy range 0.4-8.0 keV. In the first 17 fits (first 17 epochs), we fit images A, B, C and 
D with the same power-law indices but different absorption. In fits 18, 19 and 20 (corresponding to epochs 18, 19, 
and 20), we extract spectra from a circular region of radius ~3" containing images A, B, C and D, and fit them with 
a simple power-law model modified by absorption. In fit 21, we stack the last three observations, and fit the spectra 
of images A, B, C, D with the same power-law indices. In fits 22-25, we fit the stacked (20 epochs) spectra of images 
A, B, C, and D independently. In fit 26, we fit the stacked (20 epochs) spectra of images A, B, C, and D with the 
same power-law index. In fit 27, we extract and stack (20 epochs) the spectra from a circular region of radius ~3" 
containing images A, B, C and D, and fit it with a simple power-law model modified by absorption. 

b Flux is estimated in the 0.4-1.3 keV band (observed 1.3 keV corresponds to 3.5 keV in the rest frame). 

c Flux is estimated in the 1.3—8.0 keV band. 

d P(x 2 l v ) is the probability of exceeding x 2 with u degrees of freedom. 
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Table 13. X-ray Emission Line Properties. 



Lens 




Line 






Image 




Elino a 

(keV) 


(keV) 


EW a 
(keV) 


Significance 


QJ 0158-4325 




FeKa 
FeKa 
FeKa 






A 
B b 
A + B 




6 45+°'" 
6.54+g-gg 


< 0.19 

< 0.18 


64+2-11 

" u — 0.32 
0.41+0-^ 


99.5% 
99.9% 








FeKa 






A 




6-84 + oio 

— U. 1U 


< 0.12 


0.40+°^ 

— U.3U 


99.3% 




HE 1104-1805 




FeKa 






B b 


















FeKa 






A + B 




6.80+81° 

— U. 13 


< 0.29 


0.23+gl 7 


99.8% 








FeKa 




A 


+ C + D 


C 


u u — 0.17 


< 0.37 


48+2^8 

" u — 0.30 


99.5% 








FeKa 




A 


+ C + D 


C 


7 i 9 +0.07 
— U.UD 


< 0.14 


0.73+g-X 

— U.4o 


99.996% 




SDSS 0924+0219 




FeKa 






B b 


















FeKa 




A + 


C + D + 


B c 


5.81+°-f? 

— 0.15 


< 0.5 


0.28+°l* 

— 0.18 


97.9% 








FeKa 




A + 


C + D + 


B c 


7 19 +o.08 

10 


< 0.16 


o-43+g-y 


99.9% 








FeKa 






A 




6.23+1!' \i 


0.48+lif' 

^ ^— 0.17 


0.93+2^ 

v-^<->_ 0.30 


99.9994% 








FeKa 






B 




ft Kc+0.11 


< 0.15 


0.49lg;g 


99.99% 




SDSS 1004+4112 




FeKa 






C 




6.37+n f. 


0.53+°™ 

— U.zo 


0.77+2'qn 

— U.oU 


99.99% 








FeKa 






D 




6 S1+ 08 
D - ol -0.08 


< 0.12 


46+ ' 22 
U.4D_ 25 


99.89% 








FeKa 




A + 


B + C + 


D 


o - a '-0.05 


< 0.20 




1 - 3 x 10" 


-9 






FeKa 






A 




6 53+ un 
D - o,:) -0.13 


< 0.18 


U - yi -0.30 


99.3% 








FeKa 






B 




f, 79 +0.08 
' -0.08 


< 0.36 


o.9ilg;|| 


99.7% 




HE 0435-1223 




FeKa 






C 




fi S3+ 010 
"•^-O.ll 
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a Ei; ne , cTi; nc , and EW represent respectively the line position, a parameter of the Gaussian profile, and the 



equivalent line width (rest frame) with respect to the X-ray continuum. 
b The FeKa line is not detected in this image. 

c We detect a redshifted and a blucshiftcd wing of the iron line for this image. 



Table 14. Microlensing Statistics for SDSS 0924+0219 



band 


dof 


B/A 


C/A 


D/A 


B/D 


C/D 


C/B 






X 2 (prob) 


X 2 (prob) 


X 2 (prob) 


X 2 (prob) 


X 2 (prob) 


X 2 (prob) 


full 


5 


20.8 (0.0009) 


5.3 (0.38) 


4.7 (0.45) 


4.6 (0.47) 


4.9 (0.43) 


10.5 (0.06) 


soft 


5 


14.1 (0.015) 


5.8 (0.32) 


4.2 (0.47) 


6.7 (0.25) 


2.2 (0.18) 


4.5 (0.48) 


hard 


5 


7.5 (0.19) 


6.5 (0.26) 


8.7 (0.12) 


7.3 (0.20) 


5.7 (0.34) 


9.0 (0.11) 
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Table 15. Microlensing Statistics for SDSS 1004+4112 



band 


dof 


B/A 


C/A 


D/A 


B/D 


C/D 


C/B 






X 2 (prob) 


X 2 (prob) 


X 2 (prob) 


X 2 (prob) 


X 2 (prob) 


X 2 (prob) 


full 


4 


2.6 (0.38) 


19.8 (0.0005) 


33.4 (1.0 x 10- 6 ) 


144.8 (2.6 X 10- 30 ) 


44.4 (5.2 X 10~ 9 ) 


25.4 (0.00004) 


soft 


4 


5.3 (0.26) 


8.6 (0.07) 


15.7 (0.003) 


72.2 (7.7 x 10- 15 ) 


46.8 (1.7 x 10- 9 ) 


14.6 (0.005) 


hard 


4 


1.9 (0.24) 


11.1 (0.025) 


16.8 (0.002) 


73.1 (5.0 x 10- 15 ) 


11.9 (0.018) 


12.1 (0.017) 







Table 16. 


Microlensing 


Statistics for HE 0435- 


1223 




band 


dof 


B/A 
X 2 (prob) 


C/A 
X 2 (prob) 


D/A B/D 
X 2 (prob) x 2 (prob) 


C/D 
X 2 (prob) 


C/B 
X 2 (prob) 


full 

soft 
hard 


3 
3 
3 


40.8 (7.3 x 10~ 9 ) 
17.3 (0.0006) 
29.8 (1.5 x 10- 6 ) 


49.5 (1.0 x 10- 10 ) 

13.8 (0.003) 

32.9 (3.4 x 10- 7 ) 


36.4 (6.3 x 10~ 8 ) 1.4 (0.30) 
29.2 (2.1 x 10~ 6 ) 6.0 (0.11) 

13.5 (0.004) 6.9 (0.07) 


0.5 (0.08) 
3.9 (0.27) 
7.8 (0.05) 


0.7 (0.14) 
0.96 (0.19) 
0.31 (0.04) 



Table 17. FeKa Lines and X-ray Luminosities 



Lens 


Enne (keV) 


EW 


(keV) 


L (10 


44 

ergs s 


-1) a 












0.2-2 


2-10 


10-50 


QJ 0158-4325 


6.31 


±0.23 


0.55 


±0.39 


1.89 


1.51 


1.69 


SDSS 0924+0219 


6.40 


±0.40 


0.76 


±0.47 


0.80 


0.39 


0.25 


SDSS 1004+4112 


6.37 


±0.17 


0.66 


±0.27 


3.19 


2.66 


2.76 


HE 0435-1223 


6.53 


±0.10 


0.87 


±0.41 


1.68 


1.51 


1.84 


Q 2237+0305 


6.38 


±0.14 


0.60 


±0.22 


2.50 


2.37 


2.52 


HE 1104-1805 


6.90 


±0.30 


0.38 


±0.27 


4.87 


4.78 


5.96 



a X-ray luminosities in 0.2-2 keV, 2-10 keV, and 10-50 keV band (rest 
frame) integrated from the best fit spectral model corrected for foreground 
absorption. 
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Fig. 1.— Stacked images of QJ 0158-4325, SDSS 0924+0219, SDSS 1004+4112, 
HE 0435-1223, HE 1104-1805, and Q2237+0305. Note the extended emission from the 
cluster lens SDSS 1004+4112. 
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Fig. 2. — Evolution of the spectral power-law index T with time for QJ 0158—4325, 
SDSS 0924+0219, HE 0435-1223, SDSS 1004+4112, Q 2237+0305 and HE 1104-1805. For 
clarity, the first points for SDSS 0924+0219, SDSS 1004+4112 and HE 1104-1805 are shifted 
by +1500 days, and the first two events of Q 2237+0305 are shifted by +1000 days. In each 
panel, the dashed horizontal line is the best fit T to the measurements from the individual 
epochs, and the red point is the power-law index T obtained from fitting the stacked spectra. 
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Fig. 3. — Stacked spectra of QJ 0158—4325 (observed frame) fit with a simple power-law 
plus Gaussian line model modified by Galactic and lens galaxy absorption. The top, middle, 
and bottom panels show results for image A, image B, and their combined spectrum, in each 
case showing the data and the model fit in one sub-panel, and the statistical residuals in a 
second sub-panel. 
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Fig. 4. — Stacked spectra of SDSS 0924+0219 (observed frame) fit with a simple power-law 
plus Gaussian line model modified by Galactic and lens galaxy absorption. The top, middle, 
and bottom panels show results for image A+C+D, image B, and their combined spectrum, 
in each case showing the data and the model fit in one sub-panel, and the statistical residuals 
in a second sub-panel. We detect both a blueshifted and a redshifted wing of the iron line. 
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Fig. 5. — Stacked spectra of HE 1104—1805 (observed frame) fit with a simple power-law 
plus Gaussian line model modified by Galactic and lens galaxy absorption. The top, middle, 
and bottom panels show results for image A, image B, and their combined spectrum, in each 
case showing the data and the model fit in one sub-panel, and the statistical residuals in a 
second sub-panel. 



-42 - 




0.5 1 2 5 

Energy (keV) 



Fig. 6. — Stacked spectra of SDSS 1004+4112 (observed frame) fit with a simple power-law 
plus Gaussian line model modified by Galactic and lens galaxy absorption. The individual 
A, B, C, and D spectra as well as the stacked A+B+C+D spectrum are shown, in each case 
showing the data and the model fit in one sub-panel, and the statistical residuals in a second 
sub-panel. 
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Fig. 7. — Stacked spectra of HE 0435—1223 (observed frame) fit with a power-law plus 
Gaussian line model modified by Galactic and lens galaxy absorption. The individual A, 
B, C, and D spectra as well as the stacked A+B+C+D spectrum are shown, in each case 
showing the data and the model fit in one sub-panel, and the statistical residuals in a second 
sub-panel. 
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Fig. 8. — Stacked spectra of Q 2237+0305 (observed frame) fit with a simple power-law 
plus Gaussian line model modified by Galactic and lens galaxy absorption. The individual A, 
B, C, and D spectra as well as the stacked A+B+C+D spectrum are shown, in each case 
showing the data and the model fit in one sub-panel, and the statistical residuals in a second 
sub-panel. Iron line is detected in images A, B, C and D. The nickel line is only detected in 
images A and D. 
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Fig. 9. — Differential absorption n^ — n^ in the lens as a function of the power-law index 
T for QJ 0158—4325. The stacked spectra of image A and B are simultaneously fit with 
the same power-law index but different amounts of neutral absorption njj. The red solid, 
green dashed, and blue dotted curves show the two-parameter 68%, 90% and 99% confidence 
contours respectively. 
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Fig. 10. — Lens absorption nn as a function of the power-law index T for HE 1104-1805. 
The curves show the 68% confidence level contour for two parameters. 
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Fig. 11.— Absorption iih by the lens as a function of the power-law index T for 
Q 2237+0305 (upper left), HE 0435-1223 (upper right), SDSS 1004+4112 (lower left), and 
SDSS 0924+0219 (lower right). The stacked spectrum of each image was separately fit with 
a power-law model modified by Galactic and lens galaxy absorption. In SDSS 0924+0219, 
the spectra of image A, C and D were combined into one spectrum (labelled by ACD) for 
these fits. Contours are for two parameters. 
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Fig. 12.— Rest frame equivalent widths of the FeKa lines of QJ 0158-4325, Q 2237+0305, 
HE 1104-1805, HE 0435-1223, SDSS 0924+0219, and SDSS 1004+4112. Dashed lines are 
the best fits for the equivalent widths of the images as constants. For SDSS 0924+0219, we 
extracted one spectrum for the combination of images A, C, and D, and the iron line for this 
spectra is labeled by ACD (A is the brightest image). The iron line was detected in all four 
images of SDSS 1004+4112, and Q 2237+0305. We did not detect the iron line in image B 
of QJ 0158-4325, SDSS 0924+0219, or HE 1104-1805, or image D of HE 0435-1223. 
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Fig. 13. — Average rest frame FeKa line centers and equivalent widths for the six lenses. 
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Fig. 14. — Chromatic X-ray microlensing light curves for QJ 0158—4325 in the full (0.4-8.0 
keV, red circles), soft (0.4-1.3 keV, green triangles), and hard (1.3-8.0 keV, blue diamonds) 
X-ray bands. The observation dates for the soft and hard bands are slightly offset for clarity. 
Uncertainties are 1 a. 
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Fig. 15. — Chromatic X-ray microlensing light curves for HE 1104—1805 in the full (0.4-8.0 
keV, red circles), soft (0.4-1.2 keV, green triangles), and hard (1.2-8.0 keV, blue diamonds) 
X-ray bands. Uncertainties are 1 a. The observation dates for the soft and hard bands are 
slightly offset for clarity. The first observation (made in 2000) was shifted by +1,900 days 
to better resolve other observations. 
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Fig. 16. — Chromatic X-ray microlensing light curves for SDSS 0924+0219 in the full (0.4- 
8.0 keV, red circles), soft (0.4-1.2 keV, green triangles), and hard (1.2-8 keV, blue diamonds) 
X-ray bands. The first observation (made in 2005) in each panel was shifted by +1,500 days 
in order to better resolve the next 5 observations (all made in 2010). 
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Fig. 17. — Chromatic X-ray microlensing light curves for SDSS 1004+4112 in the full (0.4- 
8.0 keV, red circles), soft (0.4-1.1 keV, green triangles), and hard (1.1-8 keV, blue diamonds) 
X-ray bands. The first observation (made in 2005) in each panel was shifted by +1,000 days 
in order to better resolve the next 4 observations (all made between 2010 and 2011). 
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Fig. 18. — Chromatic X-ray microlensing light curves for HE 0435—1223 in the full (0.4-8.0 
keV, red circles), soft (0.4-1.3 keV, green triangles), and hard (1.3-8 keV, blue diamonds) 
X-ray bands. In the last row, panels (BC)/A and (BC)/D show the light curves of the flux 
ratios between the intrinsic variability template (constructed from images B and C assuming 
their variability is intrinsic) and images A and D. 
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Fig. 19. — X-ray microlensing variability amplitudes of SDSS 0924+0219 in the full, soft and 
hard X-ray bands. 
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Fig. 20. — X-ray microlensing variability amplitudes of SDSS 1004+4112 in the full, soft and 
hard X-ray bands. 
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Fig. 21. — X-ray microlensing variability amplitudes of HE 0435—1223 in the full, soft and 
hard X-ray bands. Here U (BC)" indicates the intrinsic variability template constructed from 
images B and C. 
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Fig. 22. — Equivalent width of Fe Ka lines as a function of rest frame X-ray luminosity (10-50 
keV) for SDSS 0924+0219, QJ 0158-4325, HE 0435-1223, Q 2237+0305, SDSS 1004+4112, 
and HE 1104-1805 (left to right, blue diamonds). H 1413+117 (Chartas et al. 2007) and 
MG J0414+0534 (Chartas et al. 2002) are marked by cyan squares. The black triangles (red 
circles) are AGN with absorption column densities greater (less) than 10 22 cm -2 (Fukazawa 
et al. 2011). The red line is the best fit to the low absorption systems (red points) as a 
power-law with additional intrinsic scatter in the EWs. 



